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to  the  Chief  of  the  Hydraulics  Laboratory. 
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AN  EXPERIMENT  IN  MIXING  AND  INTERFACIAL  STRESS 


PART  I :  INTRODUCTION 

1.  The  matter  of  mixing  was  one  of  the  problems  mentioned  in  the  orig¬ 
inal  letter  from  the  Chief  of  Engineers  to  the  Director  of  National  Bureau  of 
Standards,  19  December  1945,  requesting  studies  on  density  current.  Subse¬ 
quently,  a  report  on  mixing  in  arrested  saline  wedges  was  issued  (Keulegan 
1955).  The  transport  rate  of  saline  waters  into  flowing  fresh  water  above, 
Um1  ,  was  shown  to  be  proportional  to  U-Uc  ,  where  U  is  the  average  ve¬ 
locity  of  flowing  fresh  water  and  Uc  the  threshold  critical  velocity  asso¬ 
ciated  with  the  initiation  of  mixing.  Suspecting  that  the  simple  result  thus 
obtained  may  be  lacking  in  general  application,  it  was  decided  to  examine  the 
corresponding  problem  for  the  case  where  the  upper  fluid  is  flowing  over  a 
long  pool  of  mainly  stagnant  saline  water. 

2.  Although  the  experimental  part  of  this  study  was  completed  in  the 

spring  of  1960,  only  recently  the  opportunity  was  presented  to  examine  and  to 

interpret  the  observational  data.  The  main  result  is  that  the  entrainment 

ratio  U  /U  is  practically  independent  of  the  densimetric  Froude  number, 
m 

This  is  at  variance  with  the  recent  investigations  showing  that  this  ratio 
instead  of  being  a  constant  is  a  function  of  densimetric  Froude  number.  This 
disparity  possibly  can  be  related  to  the  flow  end  conditions.  In  this  case, 
the  present  study  might  be  of  value  and  accordingly  the  procedures  followed 
are  given  in  great  detail  in  order  to  facilitate  the  examination  of  the 
problem  by  other  interested  researchers. 
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PART  II:  EXPERIMENTAL  SETUP 


3.  The  lower  downstream  portion  of  a  long  flume  23  cm  wide,  connecting 
the  forebay  with  the  sea  reservoir,  was  modified  to  form  the  experimental 
setup  of  36.5-m  length.  Previously,  the  channel  was  employed  to  study  the 
properties  of  arrested  saline  wedges.  The  entrance  segment  of  the  modifica¬ 
tion  is  shown  in  Figure  1;  the  exit  segment,  in  Figure  2.  The  drawings  are 
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Figure  1.  Entrance  segment  of  experimental  channel 
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Figure  2.  Exit  segment  of  experimental  channel 

not  to  true  scale.  The  sketch  of  the  first  is  self-explanatory.  In  the 
second,  a  firm  barrier  plate  separates  the  exit  segment  from  the  sea  reservoir 
which  in  the  present  study  served  as  the  source  of  saline  waters.  Fresh  water 
leaves  the  channel  laterally  through  the  exit  orifices  located  on  one  of  the 
channel  sidewalls.  There  are  eighteen  1/2-in. -diam  (12.7  mm)  holes  that  can 
be  manipulated  to  maintain  a  water-surface  elevation  of  30.3  cm,  measured  from 
channel  bottom,  independent  of  freshwater  discharge.  In  the  manipulation, 
an  open  slit  on  one  of  the  channel  walls,  40  cm  long  and  2  cm  deep  with  its 
lower  brim  30  cm  above  the  channel  floor,  is  of  considerable  value.  Thus 
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in  all  the  runs  of  the  study  the  total  depth  of  the  liquids,  fresh  and 
saline,  at  the  segment  is  maintained  at  30.3  cm.  The  discharging  water  is 
isolated  from  the  saline  waters  below  by  means  of  the  separating  block  B 
firmly  attached  to  the  barrier  plate.  During  the  tests,  with  increasing 
freshwater  velocities  the  upstream  end  of  the  interface  would  be  depressed 
and  the  downstream  end  at  the  exit  segment  would  be  elevated.  In  the  tests 
chosen,  the  interfaces  remained  at  the  downstream  end  just  above  the  rounded 
tip  of  the  separating  block. 

4.  The  saline  layer  of  the  runs  was  prepared  in  the  following  manner. 

The  channel  was  first  filled  with  fresh  water  to  a  depth  of  30  cm  and  next, 
the  saline  water  was  introduced  laterally  through  small  apertures  at  one  of 
the  channel  walls.  The  apertures,  eight  in  number,  were  distributed  evenly 
along  the  channel  and  in  line  with  the  bottom  of  the  channel.  The  connections 
with  the  saline  reservoir  were  through  small  diameter  hoses  with  capillary 
restrictions.  The  liquid  surfaces  in  channel  and  saline  reservoir  being  kept 
at  the  same  level,  the  driving  force  for  the  flow  in  the  hoses  arises  from 
the  density  difference  of  the  two  liquids.  Thus  saline  water  is  allowed 

to  slip  slowly  into  the  lower  layer  of  fresh  water  in  the  channel.  The  flow 
is  maintained  until  the  interface  reaches  a  level  about  15  cm  high.  The  inter¬ 
face  finally  obtained  is  far  from  being  clear  and  distinct.  The  saline  water 
is  colored  and  as  the  channel  walls  formed  from  lucite  are  transparent ,  the 
condition  of  the  interface  can  be  examined  by  the  eye  with  great  ease.  Intro¬ 
ducing  a  weak  freshwater  discharge,  affecting  the  efflux  through  small 
apparatus  at  one  of  the  walls  in  the  area  of  the  tip  of  block  B,  the  interface 
is  made  sufficiently  clear  to  observe. 

5.  The  fresh  water  for  the  tests  emanates  from  a  small  constant-level 
tank  with  numerous  efflux  orifices  at  its  base,  placed  over  the  forebay  lead¬ 
ing  to  the  channel.  The  orifices  consist  of  short  brass  tubes  of  various 
sizes  and  are  calibrated  in  place  by  noting  the  time  required  to  fill  a 

channel  segment  of  a  given  length  to  a  desired  depth.  For  the  study,  eight 

3 

distinct  discharges  are  employed,  the  least  being  850  cm  /sec  and  the  largest, 

3 

7,692  cm  / sec. 

6.  To  prevent  the  excessive  lowering  of  the  interface  during  the  ver¬ 
tical  traverses,  either  of  velocity  or  of  density,  due  to  the  entrainment  of 
saline  waters,  first  the  region  under  block  B  was  connected  with  the  saline 
reservoir  through  a  short  rubber  hose  of  1/2-in.  (12.7  mm)  internal  diameter. 
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Next,  as  this  arrangement  was  not  found  to  be  effective,  the  entrance  area  of 
the  saline  waters  in  the  channel  was  connected  with  the  saline  reservoir 
through  three  rubber  hoses  of  1/2-in. -diam  tubes.  The  driving  force  for  the 
flows  through  the  hoses  is  due  to  the  differences  of  the  vertically  arranged 
densities  of  the  liquids  in  the  channel  and  the  reservoir  saline  water.  The 
flows  through  the  hoses  enter  the  channel  laterally. 

7.  A  rectangular  trough,  5  cm  wide,  3  cm  high,  and  33  m  in  length,  is 
attached  horizontally  to  one  of  the  walls  of  the  channel.  This  is  used  to 
measure  fall  of  free  surface.  It  is  connected  with  the  channel  at  midstation, 
and  the  water  in  the  trough  is  of  the  same  elevation  as  the  channel  surface 
water  at  midstation  at  all  times.  The  surface  fall  is  determined  by  means 
of  optical  manometers  to  be  discussed  later. 
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PART  III:  PROCEDURES  AND  MEASUREMENTS 


Slope  of  Interface 


8.  Interface  heights,  hg  ,  are  observed  visually  against  10  paper 
scales  attached  to  the  outside  surface  of  one  of  the  channel  walls  at  equal 
intervals.  Due  to  parallax,  the  readings  could  be  in  error  as  much  as  2  mm. 
This  is  ignored.  Initially  the  interface  is  horizontal,  but  soon  after  fresh 
water  of  a  chosen  Q  is  introduced,  the  interface  at  the  entrance  segment 
commences  to  fall.  This  is  due  to  the  fact  that  in  the  main  the  saline 

water  is  stagnant,  and  the  frictional  force  of  the  moving  water  is  subsequently 
balanced.  In  the  runs,  the  interface  height  of  the  exit  segment  tended  to 
each  a  constant  value  irrespective  of  the  freshwater  velocities.  For  saline 
water  relative  densities  Apg  =  0.02,  0.04,  and  0.06  ,  the  heights  were  17.1, 
18.9,  and  17.2  cm,  80  cm  upstream  of  the  separating  block.  Definitely,  the 
saline  height  of  the  area  is  controlled  by  the  geometry  of  the  exit  segment. 

At  the  entrance  segment,  saline  water  heights  decrease  with  increasing  fresh¬ 
water  discharges.  Neither  of  the  two  ends  is  a  critical  section. 

9.  For  small  freshwater  discharges,  saline  water  height  increases 
uniformly.  For  large  discharges,  the  increase  of  height  is  not  uniform  and 
the  rate  of  increase  is  augmented  with  distance.  Data  of  Figure  3  are  an 
example.  Whether  this  behavior  is  related  to  the  actual  flow  mechanism  or 

is  brought  about  for  difficulty  of  observation  of  a  diffused  interface  picture 
due  to  severe  turbulence  is  not  clear.  In  this  case,  the  initial  rate  of 
increase  is  taken  to  represent  the  interface  slope.  Data  on  the  interface 
slope  in  the  tests  where  both  the  fall  of  surface  waters  and  the  interface 
heights  are  observed  are  collected  in  Table  1.  The  slope  of  interface  is 
related  to  densimetric  Froude  number  as  shown  in  Figure  4.  Data  are  based  on 
evaluations  entered  in  Table  2.  Some  scatter  of  points  is  present  with  the 
data  for  Apg  =  0.04  being  displaced  higher  than  the  rest.  Yet  disparities 
are  not  systematic  and  the  straight  line  drawn  yields  the  relation: 


dh 
_ s 

dx 


U 


ADs 

- —  gh 
p  w 


C  =  1.11  x  10 


Here,  U  is  the  water  mean  velocity  computed  from 


(1) 


7 


Figure  3.  An  example  of  variation  of  interface  height  with  distance 
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u/U 


Figure  13.  An  example  of  velocity  distribution  in  saline  water  layer 
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Figure  12.  An  example  of  velocity  distribution  in  freshwater  layer 


In  the  actual  use  of  a  ribbon,  however,  ,  instead  of  being  based  on  this 

1/2 

last  relation,  was  determined  by  using  the  summation  of  6  readings  in  a 

traverse  from  the  freshwater  surface  down  to  the  interface  in  accordance  with 


h 

w 


o 


Here  z'  is  the  distance  measured  from  the  interface,  h  the  depth  of 

w 

freshwater  layer,  and  U  the  current  mean  velocity.  Small  velocities  of  the 
saline  water  layer  were  measured  mostly  using  lucite  ribbons  while  phosphor 
bronze  ribbons  were  reserved  for  the  freshwater  velocities. 

Fresh  and  Saline  Water  Velocities 

18.  Examples  of  fresh  and  saline  water  velocity  measurements  are  shown 

in  Figures  12  and  13,  respectively.  Drawing  smooth  curves  through  the  data 

points  for  a  selected  abscissa  u  may  be  read  to  form  the  ratio  u/U  ,  U 

being  the  current  mean  velocity  obtained  as  U  =  q/h  .  The  relative  veloc- 

w 

ities  of  the  freshwater  area,  from  the  runs  with  Ap  /p  =  0.04  and  observed  at 

s 

x  =  2,000  cm,  are  shown  in  Table  3.  It  is  evident  that  the  pattern  of  veloc¬ 
ities,  u/U  as  a  function  of  z'/h  ,  is  hardly  affected  by  the  current 

w 

velocity  in  the  range  from  5.5  cm/sec  to  14.10  cm/sec .  The  averages  are 

shown  in  the  last  column.  Similar  behavior  is  noted  in  the  measurements  of 

the  runs  with  Ap^/p  =  0.02  and  Apg/p  =  0.04.  The  averages  from  all  these 

tests  are  shown  in  Table  4.  There  are  only  small  differences  in  the  entries 

for  a  given  z'/h  and  thus  it  would  be  appropriate  to  take  the  averages 
w 

shown  in  the  last  column  as  the  characteristic  pattern  of  freshwater  velocities 
independent  of  the  mean  current  velocity  and  the  saline  water  density.  The 
graph  of  the  average  is  shown  in  Figure  14  and  the  latter  will  be  utilized  in 
the  entrainment  analysis. 

19.  The  values  of  u/U  in  the  saline  water  area  from  the  corresponding 
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Figure  11. 


Dependence  of  deflection  of  immersed  ribbon  on  current  veloci 


dc  d?  dc  (19) 


O  0  0  5 


since  n  and  ^  vanish  at  C  =  0  ,  due  to  the  manner  ribbon  is  attached  to 
d5 

gage  vertical  stem.  Obviously,  a  first  approximate  solution  is: 


n  = 


(20) 


and  the  analysis  may  be  continued  evaluating  the  second  and  the  third  terms 
using  the  approximate  value  of  n  .  The  rest  is  just  so  much  algebra;  the 
steps  will  be  omitted  here.  Finally,  the  end  reflection  on  this  basis  is 
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j  =  0.25II  -  0. 02511 


(21) 


17.  During  the  spring  of  1956,  H.  L.  Carpenter,  a  previous  colleague  of 
mine,  studied  experimentally  the  behavior  of  ribbons  subjected  to  currents 
using  the  same  channel  as  of  the  present  investigation.  The  ribbons  selected 
were  of  phosphor  bronze  and  Duralumin  material,  21.3  cm  in  length,  of  l/4-in. 
width  and  of  various  thicknesses.  The  ribbons  were  placed  in  currents  with 
water  depths  from  3  to  5  cm  and  at  such  distances  from  the  channel  bottom  that 
the  point  velocities  would  be  the  same  as  the  mean  velocity.  Such  distances 
were  determined  on  the  basis  of  Blasius  velocity  law.  The  scale  to  read 
deflections  was  pasted  on  the  outside  wall  of  the  channel.  Zero  of  scales 
was  read  with  water  still  in  the  channel.  Results  from  this  study  are  shown 
in  Figure  11  together  with  the  theoretical  determination  evaluated  from  Equa¬ 
tion  21.  The  alignment  of  the  experimental  points  is  linear  leading  to  a 
straight  line  slightly  above  the  theoretical  line.  The  former  may  be  de¬ 
scribed  as 
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This  suggests  that  the  relation  between 

u  =  tc  6 


6  and 

1/2 


u  may  be  expressed  as 


(23) 
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The  flexural  deformation  of  ribbon  is  in  accordance  with 


M  = 


bEt 

12R 


3 


where  E  is  the  modulus  of  elasticity  of  ribbon  material,  t  the  ribbon 
thickness,  and  R  the  radius  of  curvature  at  s  , 


1 

R 


Since  it  was  assumed  that  dy/ds  is  small,  the  latter  becomes 

I  =  tLJL  +  11  (&\ 

R  ,2  6  ds  \ds/ 

ds 

The  combination  of  Equations  12,  13,  and  15  leads  to 
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Introducing 
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Equation  16  transforms  to 


Integrating  the  terms  of  the  equation  with  respect  to  r,  twice,  between 
C  =  0  and  c  =  C 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

limits 
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Figure  10.  Deflection  of  ribbon  in  flowing  liquid 

The  total  moment  on  the  cross  section  at  s  due  to  the  moment  of  the 
forces  acting  as  the  part  of  ribbon  between  s  and  t  ,  i  the  length  of 
ribbon,  would  be 


M  =  pbCD  H_  (Xi  -  x)2 


In  terms  of  s  ,  since 


=/s(£) 


the  moment  is 


"  /I1  -  ©*] 


Restricting  attention  to  the  case  where  dy/ds  is  small  with  respect  to  unity 
..  DbCDu2L  „  2  .  f/dvi2  .  1 


u  -  s)2  -  a  -  s) 

r  (dt)  ds 
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curves  coalesce.  Since  H0  =  hs  +  hw  »  very  nearly»  the  curve  of  Figure  9 
suggests  the  relation 


d5 

dx 


0.012 


(7) 


The  Method  of  Velocity  Measurements 

15.  Fresh  and  saline  water  velocities  are  measured  by  noting  the  end 
deflection  of  a  thin  metallic  ribbon  held  normal  to  the  flow  direction.  The 
length  of  the  ribbon  is  just  slightly  less  than  the  width  of  channel. 

Attached  to  the  stem  of  a  measuring  gage,  it  may  be  raised  or  lowered  to  any 
position.  The  scale  for  the  deflections,  attached  to  a  second  gage,  may  be 
also  raised  or  lowered  to  any  position.  The  zero  of  the  ribbon  is  noted  after 
raising  the  ribbon  and  the  scale  just  slightly  above  the  water  surface.  It  is 
assumed  that  in  lowering  the  ribbon  and  the  scale  into  the  fresh  water,  on 
the  saline  water  the  scale  zero  reading  is  not  changed.  This  is  not  certain. 
Possibly  low  velocity  readings  are  open  to  error.  At  the  time  it  was  realized 
that  this  difficulty  may  be  eliminated  by  mounting  a  strain  gage  at  the  fixed 
end  of  the  ribbon.  However,  this  meant  a  new  effort  to  develop  the  new  device. 
Lacking  the  time  for  it  the  idea  was  relegated.  In  what  follows,  the  elastic 
theory  of  the  ribbon  will  be  developed  (Figure  10). 

16.  The  force  on  an  element  of  ribbon  of  length  ds  having  the  projec¬ 
tion  dx  and  dx'  normal  to  the  flow  is 

2 

AF  =  obCD  ~  dx  (8) 

where  p  is  the  density  of  liquid,  b  the  width  of  ribbon,  the  drag 

coefficient,  and  u  the  current  velocity.  The  coefficient  may  be 

taken  as  2  .  The  corresponding  bending  moment  on  a  section  at  a  point  of 
distance  s  measured  along  the  length  of  the  ribbon  would  be 

2 

AM  =  pbCD  ~  x'  dx'  (9) 

where  x'  is  the  distance  from  x  ,  the  abscissa  value  corresponding  to  s  . 
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x 2  ~  2,723  cm,  and  x^  =  3,080  cm.  A  more  efficient  use  of  the  manometers 
requires  that  each  manometer  is  provided  with  an  individual  telescope  and 
scale  firmly  placed.  In  this  case,  a  zero  of  scale  is  readily  obtained. 
Unfortunately,  only  one  telescope  was  available  at  the  laboratory  which 
necessitated  that  the  telescope  and  the  accompanying  scale  be  mounted  on  a 
table  and  the  whole  assembly  be  moved  from  one  manometer  location  to  the 
other,  keeping  the  distance  between  scale  and  channel  constant.  In  this  pro¬ 
cedure,  an  accurate  zero  cannot  be  had  and  the  difficulty  is  overcome  as'  fol¬ 
lows.  With  the  telescope  facing  a  manometer,  the  frontal  cell  is  connected 
with  the  trough  along  the  channel  and  the  hinder  cell  with  the  channel.  Scale 
reading  R^  is  noted.  Next  the  connections  are  reversed,  that  is,  the 
frontal  cell  is  connected  with  the  channel  and  the  hinder  cell  with  the  trough 
and  scale  reading  R?  is  noted.  In  this  way  the  indication  of  fall  of  water 
surface  in  the  channel,  R  ,  with  respect  to  the  water  surface  in  the  trough 
is  doubled,  that  is  =  AR  =  k  .  The  connection  of  the  trough  with 

channel  waters  is  believed  to  be  at  x  =  1,610  cm  and  thus  AR  would  be 
positive  for  the  first  manometer  and  negative  for  the  remaining  two. 

13.  The  readings  R^  and  R2  for  a  manometer  were  taken  repeatedly 
for  a  period  of  20  min.  Averaged  values  are  taken  as  representative.  The 
resulting  indications  k  obtained  during  the  runs  with  Ap g/p  =  0.06  are 
shown  in  Figure  7  as  an  example.  It  appears  as  if  the  rate  of  surface  fall 
6  is  far  from  being  uniform  and  increases  in  the  downstream  direction,  es¬ 
pecially  when  freshwater  discharge  is  large;  but  an  admissible  value  for  the 
surface  fall  is  obtained  as  follows.  Denoting  the  magnified  values  of  the 
surface  fall  with  respect  to  waters  in  the  trough  by  k^,  and  k^  cor¬ 

responding  to  the  manometer  positions  x^,  X£»  and  ,  respectively,  the 
effective  surface  fall  is: 


Determinations  made  in  this  manner  are  shown  in  the  last  column  of  Table  1. 

14.  Dependence  of  the  rate  of  all  d6/dx  on  the  Froude  number  based 
on  the  mean  velocity  of  flow  and  the  depth  of  freshwater  layer  is  indicated  in 
Figure  8.  Data  are  taken  from  Table  2.  The  curves  corresponding  to  three 
saline  waters  show  disparities  at  greater  freshwater  discharges.  On  the 
other  hand,  when  the  Froude  number  is  based  on  total  depth  the  individual 
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Figure  6.  Optical  manometer;  cylinder 
length  10  cm,  cell  length  16  cm 

manometer  mirror  by  S  and  the  scale  reading  by  R  ,  by  geometry 


Thus  the  magnification  of  displacment  6^  is  M  =  2(S/s)  . 

11.  Imperfect  meniscus  formation  along  the  cylinder  surfaces  must  be 
avoided  for  the  satisfactory  operation  of  the  manometers.  Although  this  will 
be  somewhat  less  significant  when  the  floating  cylinders  are  made  from  hydro- 
phobic  material,  for  added  precaution,  however,  small  amounts  of  aerosol  are 
added  to  the  water  of  the  manometer  cells.  The  value  of  the  magnification  M 
is  obtained  by  calibration  instead  of  relating  to  geometry  through  Equation  5. 
First,  water  levels  of  the  two  cells  are  equalized  by  opening  the  connection 
between  them  and  closing  the  connection  with  trough  and  channel.  Scale  read¬ 
ing  R,  is  noted.  Next,  the  connection  between  the  two  cells  is  closed  and 

1  3 

water  of  volume  5  cm  is  introduced  into  the  frontal  cell.  The  new  scale 

2 

reading  R£  is  noted.  The  basin  area  of  each  cell  is  160  cm  and  the  re¬ 
sulting  difference  in  heights  of  the  water  levels  in  the  cells  is  5/160.  Thus 
the  manometer  magnification  is  160(R2  -  R^)/5  .  For  the  three  manometers, 
calibration  yielded  the  values  44.9,  43.4,  and  44.0.  The  value  M  =  44.0  is 
adopted  as  applying  to  all  the  manometers. 

12.  The  manometers  were  located  at  the  stations  x^  =  300  cm, 


12 


-10*  dh/dx 


Figure  5.  Variation  of  interfacial  slope  with  densimetric 
Froude  number  (after  Loftquist  1960) 
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(2) 


w 


where  is  the  freshwater  depth  at  x  =  1,600  cm,  q  the  freshwater  dis¬ 

charge  per  unit  width  of  flume,  p  the  density  of  water,  and  p  +  Apg  the 
density  of  salt  water.  Ignoring  the  fall  of  surface  waters, 


h  =  H,  -  h 
w  1  s 


(3) 


As  mentioned  before,  ,  the  total  depth  of  the  two  layers  at  the  exit 

segment,  is  of  constant  value  30.3  cm  for  all  the  tests  of  the  study.  Loft- 
quist  (1960)  studied  the  mechanism  of  underflow  in  this  same  channel  where 
the  present  study  was  carried  out.  In  this  case,  the  interface  would  be  fall¬ 
ing  and  one  may  deduce  from  these  data  (Figure  5)  the  relation 


dh 
_ s 

dx 


U 


Ap 

— “  gh 
p  B  s 


C  =  1.35  x  10 


-2 


(4) 


where  U  is  the  mean  flow  velocity  of  saline  water  below  h^  ,  the  depth  of 
saline  underlayer.  Agreement  of  these  two  results  is  close. 


Fall  of  Surface  Waters 


10.  The  fall  of  surface  waters  is  measured  with  three  optical  manom¬ 
eters,  each  consisting  of  two  cells  that  may  be  connected  with  each  other, 
with  the  small  trough  running  along  the  length  of  channel,  or  with  the  main 
channel  itself.  Details  of  construction  are  shown  in  Figure  6.  Two  lucite 
cylinders,  the  ends  closed,  are  attached  to  a  rectangular  platform  with  a 
thin  stem  projecting  upward  on  which  a  galvanometer  mirror  is  attached. 
Distribution  of  the  masses  is  such  that  when  the  waters  of  the  cells  are 
in  the  same  horizontal  plane,  the  center  lines  of  the  two  cylinders  would 
be  also  in  the  same  level.  The  differential  of  the  levels  of  the  waters  in 
the  cells  would  lead  to  a  rotation  of  mirror  surface.  This  could  be  trans¬ 
lated  into  a  scale  reading  using  a  telescope  with  a  vertical  scale  next  to 
it.  Let  s  be  the  distance  between  the  cylinders  and  the  relative  dis¬ 

placement  of  the  waters  in  the  cell.  The  resulting  rotation  of  the  mirror 
surface  is  <5c/s.  Denoting  the  distance  of  the  telescope  scale  from  the 
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Figure  14,  Characteristic  curve  of  velocities  in  freshwater  layer 
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runs  are  shown  in  Table  5.  From  column  to  column  there  are  some  differences 
in  the  u/U  values  for  a  given  elevation  z'/hg  •  These  are  probably  due  to 
observation  difficulties,  mainly  from  the  uncertainties  of  the  scale  zero  to 
read  ribbon  deflections.  With  this  in  mind  it  is  appropriate  to  consider 
averages  as  shown  in  the  last  column.  This  average  together  with  those  from 
runs  relating  to  Apg/p  =  0.02  and  Ap^/p  =  0.04  are  shown  in  Table  6.  It 
would  be  appropriate  to  take  the  average  values  shown  in  the  last  column 
(plotted  in  Figure  15)  as  being  the  characteristic  velocity  pattern  for  the 
saline  water  layer  independent  of  current  velocity  and  saline  water  density. 


Figure  15.  Characteristic  curve  of  velocities  in  saline  water  layer 
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20.  The  above  discussions  relate  to  the  velocities  observed  at  station 


x  =  2,000  cm.  In  addition  to  these,  in  separate  tests,  the  velocities  at 
station  x  =  3,000  were  also  investigated.  The  main  results  from  the  fresh 
water  are  shown  in  Figure  16.  The  curve  drawn  is  reproduced  from  Figure  14. 
The  close  agreement  of  the  data  points  with  the  curve  suggests  that  there  is 
similarity  in  the  velocity  pattern  from  section  to  section.  The  case  for  the 
saline  water  area  is  shown  in  Figure  17  and  the  curve  likewise  is  reproduced 
from  Figure  15.  Complete  similarity,  here,  does  not  exist.  Data  indicate 
that  there  is  at  the  farther  downstream  cross  section  an  effective  flow  to 
the  right.  This  would  be  indicative  of  a  pressure  buildup  in  the  area  of 
channel  exit.  Owing  to  the  smallness  of  ribbon  deflection  and  of  possible 
zero  errors,  we  are  not  altogether  certain  of  this. 

Interfacial  and  Free  Surface  Velocities 


21.  The  nature  of  the  freshwater  velocities  is  an  important  item  in  the 
evaluation  of  entrainment  carried  out  subsequently.  It  was  concluded  in  the 
above  that  in  this  particular  area  the  velocity  patterns  are  similar  to  each 
other  at  different  cross  sections.  Interfacial  and  free  surface  velocities 
should  be  parts  of  this  picture. 

22.  The  surface  velocities  were  measured  by  placing  small  paraffin 
particles  on  the  central  part  of  water  surface  and  noting  the  travel  time 
over  selected  distances.  The  values  noted  in  the  runs  with  saline  waters  of 
relative  density  difference  Apg/p  =  0.06  are  shown  in  Table  7.  In  this, 

Ug/U  denotes  the  ratio  of  surface  velocity  to  the  average  current  velocity 
U  at  sections  of  relative  distance  x/L,  the  ratio  of  section  distance  to 
the  total  length  of  saline  pool  under  fresh  water.  Averaging  the  results 
obtained  with  varying  current  velocites  U  and  adding  the  corresponding  ones 
from  the  runs  with  Apg/p  =  0.02  and  0.04  are  shown  in  Figure  18a.  The 
line  drawn  yields  the  value  Ug/U  =  1.08  ;  but  the  corresponding  value  from 
Figure  14  is  0.98.  The  difference  in  these  two  values  may  be  attributed  to 
the  fact  that  the  paraffin  particles  measure  the  maximum  of  Ug  ,  whereas 
the  ribbon  measures  the  mean  root  square  of  Ug  across  the  entire  width  of 
channel . 

23.  For  the  observation  of  the  interface  velocities  a  globule  of  butyl 
phthalate  dissolved  in  xylene,  with  a  density  equal  to  half  that  of  saline 
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water  of  a  particular  run,  is  injected  by  means  of  an  eye  dropper  into  the 
water  a  few  centimeters  below  the  surface.  The  size  of  the  globule  may  be 
controlled  by  the  diameter  of  the  dropper  opening.  The  globule  diameters  fell 
within  the  range  of  3  to  6  mm.  After  the  descent  of  the  drop  is  completed, 
reaching  the  interface,  the  subsequent  motion  of  the  drop  from  one  locality  to 
another  is  timed.  Let  be  this  interface  velocity  observed  at  the 

point  x  and  U  be  the  mean  velocity  of  fresh  water  in  the  cross  section 
through  x.  Forming  the  ratio  U^/U  ,  observed  values  from  tests  with 
dp  /o  =  0.06  and  with  varying  freshwater  velocities  are  shown  in  Table  8,  as 
an  example.  For  a  given  x/L  ,  from  one  column  to  the  other,  values  reveal 
only  small  differences.  These  differences  may  be  ignored  and  the  average 
taken  to  represent  the  group.  These  U^/U  averages  and  the  similar  ones 
from  runs  relating  to  saline  water  relative  density  differences  Ap  /p  =  0.02 

and  0.04  are  shown  in  Figure  18.  There  are  small  differences  in  the  U^/U 
values  for  a  given  x/L  .  As  these  are  not  related  to  Ap^/p  ,  we  take  the 
curve  drawn  to  be  of  general  application.  Thus  independent  of  interface 
locality,  saline  water  density,  and  freshwater  current  velocity,  U^/U 
equals  0. 54 . 


Density  Probe 

24.  Densities  are  determined  directly  by  a  traversing  probe,  the 
exploring  end  consisting  of  two  parallel  copper  or  platinum  wires,  2  cm 

in  length,  0.8  mm  in  diameter,  and  the  spacing  between  them  0.5  cm.  A  small 
transformer  provides  the  electromatic  force,  6  volts,  and  the  alternating 
current  through  the  electrodes  is  measured  by  an  a-c  ammeter  with  shunts  for 
various  sensitivities. 

25.  The  salinity  of  an  electrolyte  is  related  to  the  resistance  between 
two  parallel  wires  immersed  in  the  solution.  If  the  spacing  between  the  wires 
is  small,  the  resistance  is  given  by 

R  =  log 
ttL  6e 

where  R  is  the  effective  resistance  in  ohms,  t  the  specific  resistance 
of  the  electrolyte  in  ohms  f-m,  L  the  length  of  the  individual  wires,  d 
the  diameter  of  the  wires,  and  s  the  spacing  of  them,  all  measured  in 


2S 

d 


1  - 


(26) 
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centimeters.  The  above  formula  which  should  apply  only  to  wires  of  great 
length  is  used  to  guide  the  design  of  the  probe.  The  probe  is  calibrated 
with  standard  solutions  of  known  densities.  The  above  equation  is  due  to 
Keulegan  (1949). 

26.  The  probe  body  consists  of  a  flat  rectangular  lucite  bar  1/4  in. 
thick  and  13  in.  long.  Attached  to  a  gage,  the  narrow  edge  faces  the  current. 
The  electrode  wires,  located  horizontally  below  the  bar,  are  normal  to  the 
bar.  The  vertical  movement  of  the  gage  allows  a  complete  traverse  from  the 
channel  bottom  to  the  water  surface.  In  the  area  of  the  interface,  the 
densities  are  determined  at  intervals  of  0.01  ft.  In  this,  the  determination 
of  the  exact  locale  of  the  interface  line,  Ap/Apg  =  0.5  ,  may  be  open  to  an 
error  of  0.005  ft. 


Effective  Interface  Layer  Thickness 


27.  Following  Loftquist  (1960),  the  effective  interface  layer  thickness 
may  be  defined  by  the  relation 


(r) 

\  /max 


(27) 


after  denoting  the  density  of  liquid  at  z  by  p  +  Ap  .  The  graphic  deter¬ 
mination  of  1^  is  illustrated  in  Figure  19.  The  evaluation  made  for  the 
three  saline  water  relative  densities,  Apg/p  =  0.02,  0.04,  and  0.06,  are  shown 
in  Table  9.  Values  change  only  slightly  with  current  velocities.  The 
averages  corresponding  to  the  three  saline  water  densities  are  0.89,  0.81,  and 

0.74  cm.  On  this  basis  1  amounts  to  0.81  cm.  Also  shown  in  this  table 

P 

is  hw/lp  >  the  ratio  of  depth  of  water  to  the  layer  thickness.  This  ratio 
tends  to  be  more  of  a  constant  value,  and  independent  of  current  velocity. 

Their  average  values  for  Apg/p  =  o.02,  0.04,  and  0.06  are  18.8,  21.1,  and 
22.1.  The  average  is  20.6  as 


{&) 

\  /  max 


z  =  0 


(28) 


and  introducing  the  average  value  of  fWl^ 


30 


Figure 


19.  Manner  of 
thickness  on 


evaluating  interface  layer  effective 
the  basis  of  densities,  1 

P 
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(£*). 


20.6  ,  z 


The  determinations  of  thickness  p  made  from  the  traverses  through  the  cross 

section  x  -  3,000  cm  are  shown  in  Table  10.  There  are  slight  increases 

with  current  velocity.  Ignoring  these,  the  average  values  are  1  =  1.06, 

P 

1.13,  and  0.89  corresponding  to  Apg/p  =  0.02,  0.04,  and  0.06,  respectively. 
Thus  the  average  value  would  be  1  =  1.02  cm.  Here,  also,  the  ratio  tends 

to  remain  constant  for  varying  current  velocity.  The  average  corresponding 
to  APg/p  =  0.02,  0.04,  and  0.06  are  14.3,  12.5,  and  16.3,  yielding  an  over¬ 
all  average  value 


=  14.4 


and  thus 


w  dAp 

Ap  dz 
s 


=  14.4  ,  z  =  0 


These  results  indicate  a  tendency  for  1  to  increase  with  distance  from 
entrance.  The  density  measurements  in  the  interfacial  region  are  sparse.  Had 
we  selected  copper  or  platinum  wires  1/64  in.  in  thickness  for  the  probe  and 
reduced  the  spacing  between  the  wires  to  2  mm,  it  would  have  been  possible  to 
observe  the  densities  of  the  interfacial  area  at  a  greater  number  of  points, 
assuring  the  certainty  of  '1  values.  For  an  underflow,  saline  water  flowing 
under  a  relatively  stagnant  fresh  water,  Loftquist  (1960)  found: 


=  13.7 


where  h  is  the  hydraulic  radius  of  saline  layer.  In  his  tests  h  /h  ,  h 
r  s  r  s 

the  depth  of  saline  waters,  equals  2.58.  This  gives 


=  35.3 


a  value  about  40  percent  larger  than  the  corresponding  value  noted  above  for 
measurements  at  x  =  2,000  cm. 
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28.  The  interfacial  layer  thickness  can  be  defined  also  in  terms  of 


velocities  by  the  relation 


03) 


The  graphic  determination  of  1  is  illustrated  in  Figure  20.  Evaluations 
made  for  the  three  saline  water  relative  densities,  Ap^/p  =  0.02,  0.04,  and 
0.06,  are  shown  in  Table  11.  The  values  change  only  slightly  with  current 
velocities  although  there  is  a  tendency  for  1  to  increase  with  current 
velocity.  The  averages  corresponding  to  three  saline  water  densities  are 
2.76,  2.48,  and  2.66  cm.  On  this  basis  1  amounts  to  2.63  cm.  Hence  thick¬ 
ness  based  on  velocities  is  twice  as  large  as  that  based  on  densities.  Also 
shown  in  this  table  is  the  ratio  h^/1  .  Independent  of  velocity,  it  tends 

to  a  constant  value.  For  Apg/p  =  0.02,  0.04,  and  0.06  the  averages  are  6.23, 
6.78,  and  6.22.  The  overall  average  is  6.41  yielding 


Figure  20.  Manner  of  evaluating  interface  layer  effective 
thickness  on  the  basis  of  velocities,  1^ 
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29.  Among  the  earlier  Investigators,  Ellison  and  Turner  (1959)  were  the 
first  to  note  that  entrainment  is  related  to  the  gross  Richardson  number 


R. 

1 


(Apg/p)hw 


(35) 


On  the  other  hand,  more  directly  influenced  by  the  mechanism  of  flow  in  the 
interfacial  area,  entrainment  is  related  to  the  local  Richardson  number 

-R.  =  -g-dp-/-z  ■  ,  z  =  0  (36) 

10  p ( 3u/dz ) 

The  findings  of  Ellison  and  Turner  would  imply  that  the  gross  and  local 
Richardson  numbers  differ  from  each  other  by  a  numerical  factor,  a  point  which 
can  be  readily  demonstrated  through  the  results  presented  above.  Squaring  the 
two  sides  of  Equation  34,  dividing  the  terms  of  the  resulting  equation  by 
Equation  29 


Ap  h  . 

s  w  _  _2  dp/dz 


(3u/dz)" 


Multiplying  both  sides  by  g/p 


gAp  h 
s  w 

ou2 


=  -2 


gdp/dz 
p  (du/ dz) ‘ 


The  numerical  factor  for  the  present  case  is  two. 


Fresh  and  Saline  Water  Densities 


(37) 


30.  The  pattern  of  the  density  distribution  in  a  cross  section  may  be 
normalized  by  dividing  Ap  by  Aps  and  s  the  distance  from  the  interface 

by  1  .  At  the  interface  Ao/Ap  equals  half  and  1  is  the  interface 

p  s  p 

effective  thickness  determined  on  the  basis  of  densities.  This  is  the  pro¬ 
cedure  used  by  Loftquist  (1960).  As  an  illustration  the  pattern  observed  at 
three  cross  sections  in  Run  27  is  shown  in  Figure  21.  Density  variations  are 
similar  to  each  other  in  all  the  sections.  The  line  drawn  is  according  to  the 
relat ion 


34 


le  of  density  disti 
cross  sections 


Ap 


ftp. 


_1 

2 


(‘ 


tanh 


(38) 


and  satisfies  the  following  three  conditions:  (a)  for  large  values  of  s/l^ 
positive,  Ap/.Va  vanishes;  (b)  for  large  values  of  s/1^  negative, 

Ap/Aps  equals  unity;  and  (c)  at  s  =  0,  (d/ds) (Ap/Apg)  =  1/1  .  It  is  no¬ 

ticed  further  that  in  the  observations  made  at  a  selected  cross  section  the 
pattern  of  densities  is  practically  independent  of  the  density  of  saline  water 
and  the  densimetrie  Froude  number  of  the  flow.  Thus  the  data  manifesting  only 
slight  variations  between  them  may  be  averaged.  The  summarized  values  for  the 
two  cross  sections  x  =  2,000  cm  and  x  =  3,000  cm  are  shown  in  Figure  22. 
The  line  drawn  is  also  in  accordance  with  Equation  38.  The  agreement  of  the 
saline  values  with  the  curve  is  quite  satisfactory.  On  the  other  hand,  for 
the  fresh  water  a  marked  disparity  from  the  curve  is  very  noticeable. 

31.  The  range  of  densimetrie  Froude  number  encountered  in  this  study 
is  0.1  to  0.5.  When  the  densimetrie  Froude  number  is  less  than  0.28, 

Ap / Ap ^  of  the  freshwater  area  above  the  interfacial  layer  nearly  vanishes. 
There  is  no  transport  of  salt  into  fresh  water.  On  the  other  hand,  when  F, 
is  greater  than  this  value  the  salt  transport  is  of  significant  amount  and 
furthermore  it  increases  with  distance  from  the  entrance.  In  a  given  section 
Ap/Ap  is  inversely  proportional  to  distance  measured  from  the  interface 
line.  Proportionality  factor  changes  with  cross  section  as  illustrated  in 
Figure  23.  It  is  convenient  toward  the  formulation  to  introduce  the  coef¬ 
ficient  B  defined  as 

?  =  hZ 
w 

The  plotting  in  Figure  24  leads  to 


and  that  in  Figure  25 


BL 

h 

w 


3  -4 

4  x  10  F, 


(40) 


w 


If)1  V 


(41) 


There  is  considerable  scattering  of  data  points  in  both  figures  due  mainlv  to 
the  difficulty  in  the  graphic  determination  of  P  . 


Partial  Entrainment  Ratio  U  ,  / U 

ml 


AO.  T'ne  transport  velocity  U  may  be  expressed  also  as 


where 


-6 

o 


The  quantity  aApgQ  ^  represents  the  mass  of  salt  transported  by  the  interface 
layer.  The  lower  limit  is  identified  as  the  depression  of  the  point  where  the 
relative  density  Ap/Ap  equals  0.99  as  measured  from  the  interfacial  line 
z'  =  0  .  Considerable  difficulty  is  experienced  in  deciding  just  what  the  ex¬ 
act  value  of  5  is.  This,  however,  is  not  a  serious  matter  since  at  this 
o 

point  u/U  is  small.  Nevertheless,  it  is  a  source  of  some  error.  As  regards 

the  upper  limit  1  /2  ,  it  will  be  remembered  that  it  changes  slightly  with  the 

P 

fresh  current  velocity  and  also  with  length  L  .  For  convenience,  however, 

1  /2  is  taken  to  equal  0.01  ft.  The  question  whether  Q  ^  varies  linearly 
with  or  is  a  constant  independent  of  L  cannot  be  answered  with  certainty. 

Now,  one  may  write 


Q  ,  u . 
si  m  i 

U  2  2  p 


(6b) 


This  expression,  with  Q  . 

evaluated 

on  the 

basis  of 

Equation  b4  and  1 

from  observed  data,  allows 

the  determination 

of  m  . 

Results  from  the  runs 

for  the 

three  saline  waters 

are  shown 

in  the 

following 

set  of  values  where 

m.,  and 

m  relate  to  observations  made  at 

the  cross 

sections  x^  =  2,000  cm 

and  x^ 

=  3,000  cm  ,  respectively. 

A  P  / 
s  p 

1  0 ,  cm 

P  z 

n*2 

1  - ,  cm 

P  J 

m  ^ 

0.02 

0.89 

1.26 

1.12 

1.A9 

0.0A 

0.80 

1.17 

1.17 

1.81 

0.06 

0.78 

1.13 

1.13 

1.79 
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(59) 


which  may  be  written  also  as 


r  2  "i1/3 

u  =  S  <FA  )  •  R  (60 

c  2  L  cj 

Utilizing  the  above  noted  set  of  values.  Equations  59  and  60  yield  separately 


Ap  J  ;■> 

K1 

k2 

0.02 

8.75 

8.75 

0.04 

9.28 

8.67 

0.06 

9.16 

8.62 

Below  densimetric  Froude  number  0.28,  /.V  vanishes  in  fresh  water  above 

the  interfacial  layer.  The  condition  of  densities  for  densimetric  Froude 
number  larger  than  0.28  was  previously  discussed. 

39.  Taking  the  values  of  U  ./!'  from  Tables  20,  21  and  22  relating 

mZ 

to  saline  waters  of  relative  densities.  A;  /p  *  0.02,  0.04,  and  0.06,  respec- 

2  2s 

tivelv.  These  are  plotted  against  10  F  in  Figure  31.  The  equation  of  the 
1 ine  drawn  is 


U 


0. 034  F^ 


(61) 


applicable  in  the  range  of  F.  from  0.28  to  0.55.  The  scatter  of  the  points 
in  the  region  inferior  to  0.28  may  be  attributed  to  the  fact  that  the  transport 
from  the  interface  layer  into  fresh  water,  as  mentioned  above,  is  exceedingly 
small  and  their  values  are  irregular  and  cannot  be  determined  accurately. 
Another  presentation  is  given  in  Figure  32,  where  the  equation  of  the  line 
drawn  is 


In  the  latter  presentation,  points  are  more  closely  placed  to  the  curve  drawn. 

Which  of  the  two  relations  for  U  _ /U  is  valid  actuallv  is  a  nuestion  that 

m2 

cannot  be  ascertained  in  the  present  study.  Tn  the  tests  made,  as  mentioned 
previously,  R  and  F.  are  nearly  related  to  each  other. 
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U,  cm /sec 


Figure  30.  Freshwater  component  of  saline  water  transport 
against  current  velocity 
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If  it  be  assumed  that  the  critical  velocity  be  independent  of  water  depths 


3 


Figure  29.  Volume  of  saline  water  entering  into 
current  from  the  interface  layer 

is  shown  in  Tables  17,  18,  and  19,  respectively.  The  quantities  hw  and  U 
appearing  in  these  tables  represent  the  average  values  of  depth  and  current 
velocities  over  the  length  L  . 

38.  In  Figure  30,  as  an  illustration,  we  have  plotted  Um2  against  U 

for  the  case  of  Ap  /p  =  0.04  ,  taking  values  from  Table  18.  Two  facts  may  be 
s 

readily  deduced.  First,  the  data  from  two  cross  sections  are  in  close  agree¬ 
ment  with  each  other.  Secondly,  Um2  is  practically  nil  below  a  critical 

velocity  U  .  The  latter  fact  signifies  that  no  saline  water  from  the 

c 

interface  layer  enters  into  the  flowing  fresh  water  below  the  critical.  One 
may  associate  with  the  latter  a  critical  densimetric  Froude  number  (F^)  or 
a  critical  Reynolds  number  (R)^  and  which  may  be  determined  by  plotting  l'm2 
against  F^  or  against  R  .  Examination  yields  the  following  set  of  values: 
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y  dx 


Imagine  that  the  interval  to  is  divided  into  equal  segments  of  a 

common  value  Ax  and  the  ordinate  values  ,  ^"’'^n-l  ’  ^n  are  note^‘ 

Granting  that  Ax  is  small,  y  varies  linearly  in  segment  Ax  . 


y  dx 


i=n 


(57) 


In  accordance  with  this,  the  integration  in  Equation  56  may  be  replaced  by 
the  summations 


Q 


s2 


U  E 


dz1  +  U 


E 


(58) 


For  the  evaluations  carried  out,  dz^  is  taken  equal  to  0.01  in  the  first 

summation  covering  the  interval  from  z1  =  0.01  to  z'  =  0.1  ft;  in  the 

second  summation,  dz^  is  taken  equal  to  0.15  ft  covering  the  interval 

z'  =  0.1  to  z'  =  h  ft.  For  a  given  z'  ,  Ap/Ap  is  read  from  curves 
w  s 

similar  to  those  shown  in  Figure  28.  The  corresponding  u/U  value  ls  read 

from  the  curve  shown  in  Figure  14  after  z ' /h^  ratio  is  formed.  The  Q 

values  thus  computed  are  shown  in  Table  16  and  these  are  plotted  in  Figure  29. 

Distribution  of  the  plotted  points  admits  a  linear  alignment,  practically, 

showing  that  'J  .  is  uniform, 
mz 


Partial  Entrainment  Ratio  U  „/U 

_ m2 

37.  The  observations  of  Ap/Ap  in  the  main  are  made  at  two  cross 

s 

sections  X2  =  2,000  cm  and  x^  =  3,000  cm.  In  every  instance  Um7  is 
established  first  evaluating  Q  in  the  manner  explained  above  and  then 
dividing  it  by  the  interfacial  length  L  corresponding  to  the  cross  section 
where  the  density  determinations  are  made.  The  entire  list  of  the  results 
covering  the  saline  water  relative  densities  Ap^/p  =  0.02  ,  0.04  ,  and  0.06 
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2 

illustration,  the  case  of  Aps/p  =  0.06  and  q  =  185.6  cm  /sec  is  selected. 
The  relative  density  values  in  the  three  cross  sections  are  shown  in 
Tables  13,  14,  and  15;  these  are  plotted  in  Figure  28.  To  avoid  the  con¬ 
fusion  of  plotted  points,  two  of  the  graphs  are  displaced  downward.  The 
relative  position  of  the  curves  indicates  that  the  transport  of  salt  increases 
with  distance  from  the  entrance. 


A  p /Aps 


Figure  28.  An  example  of  density  distribution 
in  freshwater  layer 

36.  An  integration  in  subsequent  evaluations  is  replaced  by  a  summation 
Consider  the  integral 
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Local  velocity  u  is  normalized  by  expressing  it  in  terms  of  the  mean 
velocity  of  current,  U  =  q/h^  •  The  analysis  of  entrainment  will  be  carried 
out  in  two  steps  by  considering  separately  the  transport  of  salt  in  the  inter¬ 
face  layer  and  the  fresh  water  above  the  interface  layer.  Toward  this  end 


and 


U 

m 


ml 


+  U 


m2 


ml 


Ap 

Ap 


dz' 


m2 


Ap  u 

Ap  U 
s 


dz' 


(52) 

(53) 

(54) 


where  U  ^  relates  to  salt  transport  into  the  interface  layer  and 
that  into  fresh  water  above  the  interface. 


m2 


to 


Uniformity  of  Saline  Transport  Velocity  U 


m2 


35.  The  transport  velocity  U  £  may  be  expressed  also  as 


m2 


^s2 

L 


(55) 


where 


<s2 


^  dz' 

Ap  U 
s 


(56) 


The  quantity  acdpsQg2  represents  the  mass  of  salt  transported  per  second 
by  the  part  of  fresh  water  lying  above  the  interface  layer  upper  limit 
z  =  1  /2  .  If  one  finds  that  Qg2  varies  linearly  with  L  ,  then  this 
will  indicate  that  entrainment  velocity  is  uniform.  The  matter  can 

be  examined  for  larger  freshwater  discharge  considering  the  Ap/Apg  values 
for  the  three  cross  sections  at  x  =  1,000,  2,000,  and  3,000  cm.  As  an 
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increase  in  the  value  of  ot  when  the  solution  is  dilute.  For  the  present, 

we  will  take  a  to  have  the  constant  value  1.60  also  in  dilute  solutions. 

34.  Density  of  the  saline  water  is  p  +  Ap  and  V  is  the  velocity 

s  m  J 

of  entrainment,  then  the  mass  of  salt  transported  across  the  interface  of 
length  L  would  be 

L 

aAp  U  dx 
s  m 


Again,  on  the  same  basis,  the  mass  of  salt  transported  by  the  fresh  water 
across  the  cross  section  at  x  -  L  would  be 

h 


aApu  dz' 


o 


The  distance  z'  is  measured  from  the  interface  line  where  Ap/Apg  =  1/2  . 

The  velocity  and  density  of  fluid  passing  the  point  z'  are  u  and 

p  +  Ap  ,  respectively.  The  upper  limit  h  represents  the  elevation  of 

w 

the  free  surface  from  the  interface  line  and  the  lower  limit  6^  represents 
the  depression  of  the  point  where  4p/Aps  equals  0.9 9.  The  determination 
of  6q  is  open  to  error,  and  it  may  be  inferred  from  Figure  22  that  is 

of  the  order 


6  = 
o 


21 

P 


(49) 


Since  water  initially  traversing  the  cross  section  at  x 
salt,  the  conservation  of  mass  requires  that 


I, 


aAp  U  dx 
s  m 


aApu  dz' 


0  is  devoid  of 


(50) 


Assuming  that  is  uniform  along  the  interface,  a  is  of  constant  value, 

the  above  yields 


U 

m 


Ap 

Aps 


(51) 
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PART  IV:  SALINE  ENTRAINMENT  INTO  FRESH  WATER 


Formula  of  Entrainment 

33.  The  transport  of  saline  water  through  the  interface  will  be  deter¬ 
mined  on  the  basis  of  the  transport  of  salt  by  fresh  water  in  a  cross  section. 

Let  a  solute  of  mass  M  and  of  density  P  be  dissolved  in  a  pure  water  of 

m 

volume  d  .  Let  Si  _  and  p  +  Ap  denote  the  volume  and  the  density,  respec¬ 
tively,  of  the  resulting  solution,  p  being  the  density  of  water.  Denoting 

the  volume  of  the  solute  by  ; l  ,  the  mass  of  solute  is  p  Q  .  As  the  mass 

m  mm 

of  the  solution  is  equal  to  the  sum  of  the  masses  of  the  components 


(p  +  Ap)Q  =  p  U  +  pfi 

s  mm 


The  corresponding  relation  for  the  volume  would  be 


a  =  D  +  esi 

s  m 


where  i  is  a  factor  to  be  determined.  From  Equations  44  and  45  we  obtain 


p  42  =  aApfl 

mm  s 


where 


a  =  (l  -  e  — 

V  Pm 


The  relation  of  a  to  salinity  S  would  be 


Usually  S  is  defined  as  a  ratio  expressing  the  number  of  grams  of  salt  per 
kilogram  of  saline  solution.  In  Table  12,  the  relation  between  a  and 
An/o  is  shown.  For  larger  salinities  a  may  be  taken  as  1.51.  In  more 
dilute  solutions  a  increases  according  to  this  tabulation.  Meanwhile, 
for  S  =  0.01  ,  the  corresponding  value  of  e  is  1.032,  a  number  greater  than 
unity.  This  is  not  very  likely  to  happen  and  therefore  we  question  the 
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to- 50L/h 


Figure  27. 


Dependence  of  coefficient  of 

on  a  new  number  F.R 

A 


B 
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to^S 


Figure  26.  Relationship  between  Reynolds  number  and 
densimetric  Froude  number  in  tests 


40 


32.  Here  we  encounter  a  difficult  situation  as  BL/h^  may  be  expressed 

as  a  function  either  of  densimetric  Froude  number  or  of  Reynolds  number.  This 

may  be  brought  about  if  because  of  limited  range  of  quantities  employed  during 

the  tests  there  is  a  definite  relationship  between  the  Reynolds  number  and 

2 

the  densimetric  Froude  number  encountered.  In  Figure  26,  the  F^  and  R 
values  noted  during  the  measurements  at  the  cross  section  x  =  2,000  cm  are 
shown.  Equations  39  and  40  imply  that 

102F2  =  0.82  (10'3R)  (42) 

and  the  curve  drawn  through  the  larger  F  values  is  in  accordance  with  this 
relation.  There  is  also  a  third  possibility  of  representation.  Again,  on  the 
basis  of  Equation  40  and  Equation  41  we  also  have 

jp  =  4.88  108F~2R-1  (43) 

W 

2 

The  observed  values  of  SL/h^  are  plotted  against  F„R  in  Figure  27.  In 
the  latter  figure,  the  dispersion  of  data  points  is  less  severe  than  those  in 
Figures  25  and  26.  Hence,  Equation  43  is  a  more  valid  relation  for  Bl./h^  . 
Apparently  it  will  be  difficult  in  density  current  research  to  evaluate 
exactly  the  role  of  densimetric  Froude  number  or  Reynolds  number  on  a  partic¬ 
ular  item  if  the  magnitude  of  the  physical  quantities  involved  fall  in  a  small 
range. 
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Figure  24.  Dependence  of  coefficient  R  Figure  25.  Dependence  of  coefficient 

on  densimetric  Froude  number  on  Reynolds  number 


densimetric  Froude  number.  The  line  drawn  in  the  figure  yields  the  relation 


=  1.5  *  10  4  (1  +  1.1  F^)  (67) 

Entrainment  Ratio  U  /U 
m 


41.  The  entrainment  ratio  U  / U  is  the  sum  of  the  partial  entrainment 

m  r 

ratios  U^/U  and  U^/U  "  The  computed  values  from  the  observations  are 

shown  in  Tables  20,  21,  and  22  for  the  saline  water  densities  Ap  /p  =  0.02  , 

2  s 

0.04,  and  0.06.  These  are  next  plotted  versus  F,  in  Figure  34.  Although 


Figure  34.  Dependence  of  entrainment  ratio  upon 
densimetric  Froude  number 
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the  distribution  of  points  suggests  a  linear  alignment,  preferring  a  power 
formula,  a  curve  is  drawn  in  accordance  with 


U 

m 

U 


3  x  10 


Suga  (j.975)  gives  the  formula 


=  6 


10  3  F^3 

A 


(68) 


(69) 


deduced  from  experiments  on  arrested  wedges  carried  out  in  a  channel  of 
100-m  length  and  80-cm  width  and  with  saline  waters  in  the  range  of 
Ap  /p  =  0.005  to  0.03  .  Here  U  denotes  the  river  velocity  at  the  tip  of 

wedge  and  is  based  on  the  velocity  of  river  and  water  depth  at  the 

tip  of  wedge.  The  range  of  river  depth  is  40  to  80  cm  and  the  discharge,  q  , 

2  2 
10  to  600  cm  /sec.  Both  formulas  give  like  values  of  U  /U  at  F,  =  15  ; 

2  mi 

at  F  =  5  ,  Suga  formula  gives  a  value  four  times  smaller  and  at 

2  a 

F,  =30  ,  four  times  greater  than  the  ones  from  Equation  67. 

42.  The  disparity  between  the  evaluation  from  Equations  68  and  69 

is  a  matter  of  much  puzzlement.  It  will  be  recalled  that  the  results  of  the 
present  study  on  entrainment  are  based  on  the  density  traverses  at  two 
cross  sections  x ^  =  2,000  and  x^  =  3,000  cm,  assuming  that  velocity  of  en¬ 
trainment  for  a  given  water  discharge  is  uniform  and  constant  across  the 

entire  length  of  the  interface.  Ellison  and  Turner  (1959)  referring  to 
test  conditions  relating  to  surface  jets  over  saline  waters  state:  "Mixing 
began  as  before,  but  the  rate  of  increase  of  depth  soon  became  smaller  and 
at  a  certain  distance  downstream,  depending  on  the  density  difference  and  the 
rate  of  flow,  the  'jet'  region  changed  smoothly  to  one  where  depth  changed 
little  with  distance.  The  turbulence  in  the  mixed  layer  was  damped  and  there 
was  no  further  mixing  with  the  salt  solution  below;  the  layer  then  remained  of 
constant  depth  until  it  reached  the  second  weir."  In  this  study,  signifi¬ 
cantly  the  interface  layer  gradually  became  more  turbulent  with  distance, 
noted  during  the  velocity  measurement  of  the  interface  by  globules. 

43.  In  the  Ellison  and  Turner  study,  mixing  between  saline  and  fresh 
waters  is  confined  to  a  small  region  next  to  flow  entrance.  Anwar  and 
Weller  (1981)  give  additional  information  and  in  particular  regarding  the 
dimensions  of  the  area.  Freshwater  depth  increases  linearly  with  distance 
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from  the  entrance.  Denoting  the  depth  of  fresh  water  at  the  entrance  by  h 

o 

and  the  length  of  the  mixing  zone  by  x  ,  the  ratio  x/h  increases  with 

o 

the  densimetric  Froude  number  F„  ,  x/h  <  3  when  F.  =  0.7  ,  x/h  =  10 

A  o  A  o 

when  F ,  =  1.85  .  When  F  equals  0.52  the  entrainment  zone  almost  disappears. 
The  maximum  F,  met  with  the  present  study  is  always  less  than  0.5,  and  if 
the  Anwar  and  Weller  findings  remain  applicable,  then  mixing  zones  at  the  flow 
entrances  should  be  absent.  The  measurement  of  the  saline  water  heights  h 
in  this  present  study  commences  with  the  cross  section  at  x  =  133  cm.  'Down¬ 
stream  of  this  cross  section,  the  depth  of  fresh  water  h  continually  in- 

w 

creases.  The  condition  of  the  interface  upstream  of  this  initial  cross 
section  is  not  studied;  it  may  be  assumed  that  in  this  area  the  interface 
remains  practically  horizontal  such  as  shown  in  Figure  35. 

44.  Anwar  and  Weller  using  data  from  various  investigations  in  addition 
to  their  findings  give  the  form 

U 

~  =  4.5  x  i(f J  (70) 

for  large  values  of  F  .  This  is  sufficiently  close  to  the  form  proposed  by 

Suga.  Results  from  Loftquist  (1960)  have  been  considered.  The  I.oftquist 

results,  however,  relate  to  an  underflow  and  the  experiments  were  carried  out 

in  the  same  channel  as  the  present  investigation.  An  important  point  to  note 

is  that  Loftquist  in  arriving  at  a  value  of  V  'V  takes  U  to  be  uniform 

mm 

all  along  the  entire  interface.  Then  it  is  very  surprising  to  find  that 

U  /U  values  would  agree  with  each  other  whether  the  determinations  are 
m 

made  from  a  short  area  of  the  interface  or  from  the  entire  interfacial  area. 

45.  Probably  the  deviation  of  U^/U  values  of  the  present  investiga¬ 
tion  from  those  of  other  investigations  is  related  to  the  end  conditions  of 
the  experimental  environment  channel  geometry  of  the  entrance  segment  as 
shown  in  Figure  1.  If  Huppert  and  Britter  (1982)  criteria  for  flow  separation 
are  applied  to  the  present  test  conditions,  since  the  flows  are  subcritical  on 
the  two  sides  of  the  inclined  entrance,  a  separation  of  flow  over  saline  waters 
is  expected.  The  resulting  flow  pattern  could  not  be  conducive  to  mixing. 

Note  also  the  unusual  manner  in  which  the  fresh  water  is  removed  from  the  exit 
segment  shown  in  Figure  2.  Theoretically,  this  matter  is  important  and  it  may 
be  desirable  to  reexamine  mixing  characteristics  of  a  channel  where  the  two 
ends  have  inclined  bottoms  and  the  flow  out  is  over  a  weir. 
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Figure  35.  Saline  water  interface  elevation  of  channel  entrance  segment 


PART  V:  INTERFACIAL  AND  BOTTOM  STRESSES 


Hydrodynamics  of  Freshwater  Layer 


46.  Let  the  x-axis  be  drawn  in  the  direction  of  the  freshwater  flow, 
z-axis  vertically  upward,  and  y-axis  laterally  with  the  limits  -b  and  b  . 
Let  B  ,  H  ,  h  ,  and  h  denote  the  width  of  channel,  the  liquid  total 
depth,  depth  of  fresh  water,  and  depth  of  saline  water,  respectively.  The 
velocity  components  along  x  ,  y  ,  and  z  are  u  ,  v  ,  and  w  ,  respec¬ 
tively  (Figure  36).  Ignoring  the  turbulent  stresses  -p  u'u'  and  -p  w'w' 


Z 


in  comparing  with  the  pressure  p  and  utilizing  the  condition  of  continuity, 
neglecting  v  , 


3u  3w 
lx  3z 


0 


(71) 


and  the  equations  of  motion  in  terms  of  Lamb's  notation  are 


3u2  3 

3x  3z 


(uw) 


_  1  AE  +  1  '  pyx  +  1  3pzx 

p  3x  p  3y  p  3z 


2  3 
u  — 


-g 


I  l£  +  I  3P*Z  + 

p  3z  p  3x 


P  3y 


(72) 

(73) 


P  being  the  density  of  fresh  water.  Here,  p  denotes  the  tangential  stress 

yx 

that  the  liquid  of  increasing  y  exerts  on  a  plane  surface  normal  to  y-axis 
and  in  the  x  direction.  Similarly,  is  the  tangential  stress  which  the 
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liquid  of  increasing  z  exerts  on  a  plane  surface  normal  to  z-axis  and  in 
the  direction  of  x  .  Analogous  definitions  apply  to  p^  and  p  recall¬ 
ing  that  the  first  letter  of  the  subscript  refers  to  the  plane  acted  upon  and 
the  second  to  the  direction  of  the  stress.  Specializing  for  the  problem  at 
hand,  in  the  absence  of  wind 


p  =0,  z  =  H  =  h  +  h 
zx  w  s 


p  =  x  ,  z  =  h 
xz  s  s 


where  x  is  the  traction  of  the  fresh  water  on  the  saline  water  below.  It 
s 

is  a  positive  quantity  and  will  be  referred  to  as  the  interfacial  stress.  De¬ 
noting  the  resistive  force  of  the  vertical  wall  by  ,  a  positive  quantity. 


x  ,  y  =  b 
w 


x  ,  y  =  -b 
w  J 


47.  Equation  73  may  be  simplified.  Multiply  the  terms  of  this  equation 
2 

by  h^/U  ,  U  is  the  freshwater  mean  velocity;  the  result  is 
2  h  h  /  3p  3p  \ 

^h  +  I  -S  _|E  +  ^JEE  +  ^£  (76) 

^2  w  3x  Vu/  °  ^2  p  ^2  \  3z  3x  3y  J 

When  the  entrainment  velocity  U  and  the  interfacial  slope  each  are  small 

,m 

quantities,  the  expression  h^  —  J  would  practically  equate  the  ratio 
of  the  freshwater  depth  to  the  radius  of  curvature  of  the  interfacial  line. 

In  channels  of  horizontal  bottom,  this  ratio  is  a  small  quantity  and  hence  the 
left-hand  term  of  the  above  equation  can  be  neglected.  Next,  introducing 
the  fluctuations 


and  writing 


v  _  w_ 

C1  ’  U  C2  ’  U  f:3 


X  z 

C  =  L  *  n  =  h 


here  L  is  the  total  length  of  the  saline  water  layer,  we  may  place 
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.  h  3p  h  _  .  h  3p 

1  w  xz  w  a  — -  .  1  w  *yz  d  — - 

p  y2  V*  =  L  ac  elc3  a!,d  P  ~2  3y  =  3c  E2E3 

where  ,  e2  ,  and  are  known  to  be  small  quantities.  Hence,  for  the 

case  glWlT  >  1  the  last  two  terms  of  Equation  76  can  also  be  neglected, 
leaving 

0  =  -r,  _  1  IP- 


which  indicates  that  Lhe  pressures  are  hydrostatic.  Integrating  with  respect 
to  7.  ,  assuming  that  the  atmospheric  pressure  is  reduced  to  zero  and  assuming 
that  entrainment  hardly  modifies  the  picture  of  the  densities,  this  yields 

p  =  pg(H  -  z) ,  z  >  hg  (78) 

The  corresponding  relation  for  the  saline  water  layer  is 

p  =  pg(H  -  z)  +  Apgg(hg  -  z)  (79) 

where  Apg  represents  the  excess  of  saline  water  density  over  that  of  the 
fresh  water. 

48.  In  the  case  that  mixing  across  the  interface  is  appreciable,  the 
pressure  term  for  the  fresh  layer  need  be  expressed  in  greater  detail.  For 
simplicity  it  will  be  assumed  that  the  densities  in  cross  sections  are 
averaged  and  the  average  value  is  p  +  Ap  ’  ,  z  >  h^  ,  and  Ap  '  is  inferior 
to  Ap  .  One  has  instead  of  the  form  in  Equation  78 


yield ing 


p  =  (p  +  Ap ' )g(H  -  z) ,  2  >  h 


1  up  3H  .  \  3  Ap 

-  r* -  =  g  r—  +  g(H  -  z)  t - 

p  3x  3x  3x  p 


Integrating  with  respect  to  z  between  limits  h  and  H 


3  p  ,  3H  ,  £  ,  2  3  Ap’ 

-■r*-  dz  =  g  h  +  o  h  r  - 

3x  3x  w  2  w  dx  p 


But  by  continuity  of  mass,  nearly 
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X 

I 


Ap  U  dx  =  UAp’h  ,  U  =  ±  U  ,  +  U  0 
s  mo  w  mo  2  ml  m2 


and  thus 


n 

I 


13p  3H  hw  ^Ps  ^mo 

p^dz  =  8^hw  +  gr  —  ~r 


(80) 


h  Ap  U 

Hence,  if  mixing  is  appreciable  the  term  -g  —j  — - jj—  need  be  introduced 

into  the  right-hand  side  of  Equation  89. 

49.  Integrating  Equation  72  with  respect  to  z  ,  between  the  limits 
H  and  hg  ,  taking  p  from  Equation  78  and  using  the  relation 

H  =  H  -  6  .  (81) 

o 


where  H  is  the  initial  value  of  H  at  x  =  0  and  6  is  the  fall  of  the 
o 

water  surface,  one  now  has 


[  jtdZ+  (UW)H 

n 


T 

S 

P 


(82) 


In  the  presence  of  entrainment 

3h 

w  =  w,  =  u .  t  — ■  +  U  (83) 

l  h  i  3x  m 

s 

where  u,  and  U  are  the  interface  and  entrainment  velocities,  respectively, 
l  m 

Also, 


w 


H 


3H 

3x 


(84) 


u  being  the  velocity  at  the  water  free  surface.  The  left-hand  member  of 
Equation  82  in  view  of  these  now  changes  to 


dz  +  u 


2  3H 
s  3x 


2 


u  . 
l 


3h 
_ s 

3x 


u.U 


l  m 


(85) 
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Introducing  the  vertically  averaged  velocity  u  at  a  given  y  and  the 


Boussinesq  coefficient  u 


_2 

a  u  h 

w 


2 

u 


dz 


h  = 
w 


h 

s 


(86) 


Differentiating  with  respect  to  x  in  accordance  with  Leibnitz  rule  and 
assuming  a  to  be  independent  of  x 


(87) 


Finally,  utilizing  Equations  85  and  87,  Equation  82  changes  to 


a 


u.U 


i  m 


g 


d6  n 
dx  w 


3p 


3y 


dz 


T 

S 

P 


(88) 


Next,  multiplying  the  terms  in  this  equation  by  dy  ,  and  integrating 
between  the  limits  b  and  -b  and  then  dividing  by  B  ,  the  channel  width, 


U.U  =  g  ~  h 
r  m  dx  w 


2t  h 
w  w 


(89) 


when  U  represents  the  freshwater  mean  velocity  over  a  cross  section  and 

U.  and  U  are  the  average  values  of  u.  and  u  across  the  channel  width, 
l  m  &  l  m 

It  is  tacitly  assumed  that  the  individual  values  u.  and  u  differ  but 

l  m 

little  from  their  average  values.  This  is  done  for  simplicity  for  otherwise 
the  aim  of  accuracy  would  require  coefficients  difficult  to  obtain.  The  total 
wall  shear  averaged  over  the  depth  h  is  denoted  by  .  In  view  of  the 

mass  flow  condition  of  continuity.  Equation  89  may  be  modified  further.  Sub¬ 
ject  to  Boussinesq  approximation 


Uh  =  q  +  U  x  (90) 

w  m 

where  q  is  the  initial  water  discharge  at  x  =  0  computed  per  unit  width 
of  channel.  Assuming  that  U  is  Independent  of  x  ,  it  may  be  shown  that 
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Since  the  slope  of  surface  water  is  negligible  in  comparison  with  the  slope 
of  interface 

3h  3h 

w  _  _  s 

3x  3x 


then  the5  last  relation  may  be  written  as 

~  (u2h  )  =  2U  U  +  U2  — — 

3x  \  w y  m  dx 

Introducing  this  in  Equation  89  and  afterward  dividing  the  terms  of  the 

2 

resulting  equation  by  U  ,  one  obtains  for  the  freshwater  layer  the  final 
form  of  the  momentum  equation 


with 


8h, 

U 


w  dS 
2  dx 


X  h 
o  w 

B 


a 


dh 
_ s 

dx 


U.\  U 

l  \  m 
U  /  U 


(92) 


X 

o 


2r 

w 


Needless  to  state,  the  last  term  in  Equation  92  shows  the  effect  of  entrainment 
on  the  momentum  equation. 

Hydrodynamics  of  Saltwater  Layer 

50.  During  the  experiments  in  order  to  maintain  a  constant  disposition 
of  the  interface,  saline  waters  were  introduced  laterally  at  the  two  ends  of 
the  channel.  This  was  meant  to  balance  the  entrainment  inflow  into  the  flow¬ 
ing  fresh  water.  The  same  result  could  have  been  effected  equally  well  if 

the  inflow  was  introduced  laterally  at  the  level  of  channel  bottom  and  equally 
distributed  along  the  entire  channel  length.  The  hydrodynamics  of  the  saline 
layer  will  be  developed  for  this  latter  case. 

51.  Again  multiplying  the  terms  in  Equation  72  by  dz  ,  but  now 
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integrating  between  the  limits  z  =  0  and  z  =  li^  and  taking  p  from  Equa- 
t ion  80  and  since  U  =  0  at  z  =  0  , 


where  i ^  is  the  bottom  stress.  This  is  a  positive  quantity  as  the  flow 
adjacent  to  bottom  is  in  the  opposite  direction  of  the  freshwater  flow.  We 
can  write  the  left-hand  member  in  the  form 


,2  ,  dh 

■Ju  ,  ,  2  s 

•  - —  d  z  +  u  .  —  +  u  .  u 

Jx  l  dx  i  m 


(94) 


As  the  average  value  of  u  in  a  cross  section  vanishes  we  resort  to  U  to 
form  an  expression  similar  to  Equation  86.  Introducing  a  new  coefficient 

h 

s 

u“  dz  ,  z  <  hg  (95) 

o 


Differentiating  with  respect  to  x 


(96) 


Utilizing  Equations  95  and  96,  Equation  93  may  be  written  as 


+  M .  U  = 
l  m 


(97) 


Multiplying  the  terms  of  this  equation  by  dy  ,  integrating  between  the  limits 

2 

b  and  -b  and  dividing  the  terms  of  the  resulting  equation  bv  U  B 
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I’ll 


d  ■ 

if  X 


dh  \ 

V 

u:  dh 

1)  r  i: 

s 

\  +  . 

■  -  -  -S  4- 

s  +  ,T1 

tlx 

) 

i:2  dx 

1  B 

..L" 


which  is  tin1  n*>nc*ntiim  equation  relating  to  tlu*  saline  lavor .  Subt  ra. 
liquation  9J  from  Equal'  ion  98  there  results 


U2  i: 


8lt 

till 

dh 

•  S  '  S  + 

•  - 

dx 

1  dx 

'•  -  +  1 

o  b  ib  r 


i 

tin1  expression  to  evaluate  the  channel  bottom  stress. 

Experimental  Evaluation  of  Stresses 

r)2.  The  evaluation  of  the  interfaeial  sfess  eoet’fii  ieut  will 
effected  nep,  I  ee  t  i  in;  the  entrainment  terra;  that  i  ,,  usinp 


s  w  d  ' 


dy.  o  I! 


dii 

dx 


,  ii“  r" 

In  the  ease  of  the  bottom  stress,  the  entrainment  tern  will  he  letup 
the  quant  i  t  v  '•  jh  /B  will  be  ignored;  that  is,  us  in}; 


vll  d  : 


qii 


V 


i  ih 


dh 

s  s 

dx 


'  w  + 


B  Id 


ih  r 

t  I 


t  is  su  f  f  1  e  i  ell  t  to  place  1'  />' 


ei.ii i a  1  to  1  .  )  10  .  Saline  liquid 

ni 

coat  act  with  tin:  vert  ieal  wall-.  a>  the  upper  and  1  owe  r  parts  rioves  i- 
dire,  lions  a:nl  this  funt  makes  ,  a  uenlicible  quant  i t  v .  ! he  ohset 

l 

■■  ■ '  the  quantities  u,  h  ,  n  ,  *  l : .  dx  ,  d  /dx  are  viven  in  !nKle 

the  basis  ot  the  data  in  ladle,  h  and  /  ,  •  nrove-,  t-i  be  I. 'll  and 

Sin,  ■  E  f ) .  1  \  ,  .  Cl,  ,1  )  equals  1.1 .  1 1 1  i  .  In  ,  on  f  o  rm  i  I  v  with  ’  ■ 

I  I  m 

,-l  111. i. in-.,  tin1  wall  stress  due  t  o  i  hi'  f  I  ..  ■  •  1  fresh  water  i  ■ 


(i 


(981 

(99  ) 

t  inn 

( inol 


(  101  ) 

•d  hut 

(  I '  ’  1 

i  n 

d  v  tint 


fresh  water 


TW  -2  /  ub  \  1/4 


,  a  =  0.0225 


where  u  is  the  average  velocity  at  z  ,  averaged  over  the  width  from 
y  =  -b  to  y  =  b  .  On  this  basis 


where 


I  (s)' 


dr,  ,  n 


Using  the  data  in  Table  4,  integration  yields  A  =  0.958  and  that 


A°  =  pu2 


=  0.0432 


Determinations  of  X  for  the  tests  are  completed  in  Table  24.  In  the  same 
o 

table  are  shown  also  X  h  /B  and  dh  /dx  .  These  quantities  are  of  the  same 

o  w  s 

order  and  as  it  cannot  be  guaranteed  that  the  friction  of  the  wall  is  deter¬ 
mined  with  the  same  accuracy  as  the  interfacial  slope,  it  would  be  preferable 


to  conduct  tests  with  wider  channels. 
53.  The  computed  values  of  i 


on  the  basis  of  Equation  31  art 


entered  into  the  third  column  of  Table  25  and  these  are  shown  in  Figure  37. 
Pedersen  (1980)  using  the  results  of  various  investigations  arrived  at  the 
following  expression  for  the  interfacial  stress  coefficient.  Writing 


f  J2  =  i  ,/p  (U  -  U.) 

i  l  m  l 


(U  -  U.)  (y  -  v  )/\ 
m  i  '  o 


the  expression  is 


2/ f  .  =  2.45 

l 


i  (r  >J7~/ 

n  \  l*  l 


2  -  1.  J 


This  relation  is  meant  to  apply  equally  to  saline  underflow  and  freshwater 
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Table  7 


Water-Surface 


Veloc it ies 


An  /p  =0.06 
s 


i 


Run 

25 

26 

27 

_ 28 _ 

29 

30 

Mean 

U ,  cm/ sec 

5.46 

3.01 

10.25 

11.76 

13.50 

15.09 

15.73 

x/L 

11 

s 

/U 

0.095 

1.144 

1 . 196 

1.206 

0.949 

1.075 

0.942 

1.030 

1.077 

0.319 

1.137 

1.196 

1.161 

1.066 

1.023 

1.090 

1.127 

1.114 

0.540 

1.067 

1.160 

1.140 

1.121 

1.094 

1.069 

1.102 

1.107 

0.768 

1.133 

1.136 

1.070 

1.059 

1.168 

1.102 

1.078 

1.099 

0.940 

1.050 

1 .048 

1.182 

1.160 

1.123 

1.127 

1.111 

1.114 

Table  8 

Interfacial  Velocities 
Ap  /o  =  0.06 


Run 

24 

25 

26 

27 

28 

29 

30 

Mean 

cm/ sec 

5.46 

8.01 

10.25 

11.76 

13.50 

15.09 

15.73 

x/L 

!! 

./u 

1 

0.095 

0.493 

0.523 

0.538 

0.  578 

0.561 

0.  570 

0.535 

0.542 

0.319 

0.552 

0.576 

0.540 

0.  566 

0.560 

0.547 

0.53  5 

0 .552 

0.540 

0.495 

0.555 

0.517 

0.515 

0.540 

0.526 

0.511 

0.507 

0.7  68 

0.4-7 

0. 524 

0.514 

0.520 

0.493 

0.533 

0.483 

0.  507 

0 . 440 

0.492 

0.501 

0.47  5 

0.516 

0.440 

0.503 

0.4  64 

0.484 

Table  6 

Velocity  Distribution  in  Saline  Water  Pool 


Ao  /p 

s 

0.02 

0.04 

u/U 

0.06 

Mean 

-z  *  /h 

s 

0.00 

0.  540 

0.540 

0.  540 

0.540 

0.02 

0.390 

0.415 

0.390 

0.  347 

0.04 

0.306 

0.286 

0.266 

0.286 

0.06 

0.258 

0.204 

0.191 

0.216 

0. 10 

0.197 

0.135 

0.124 

0.152 

0.15 

0.145 

0.079 

0.080 

0.101 

0.20 

0.108 

0.045 

0.055 

0.069 

0.30 

0.064 

0.000 

0.015 

0.026 

0.40 

0.035 

-0.030 

-0.014 

-0.003 

0.50 

-0.033 

-0.052 

-0.055 

-0.047 

0.60 

-0.043 

-0.070 

-0.061 

-0.057 

0.70 

-0.118 

-0.080 

-0.100 

-0.100 

0.80 

-0.125 

-0.079 

-0.095 

-0.100 

0.85 

-0.115 

-0.063 

-0.070 

-0.082 

0.90 

-0.098 

-0.061 

-0.068 

-0.076 

0.95 

-0.062 

-0.041 

-0.040 

-0.047 

1.00 

-0.000 

-0.000 

-0.000 

-0.000 

Table  5 


Velocity  Distribution 

in  Saline 

Water  Pool 

Ap,/p  = 

0.040; 

+  h  =  30 
w 

.  3  cm 

Run 

18 

19 

20 

21 

22 

23 

Mean 

h  ,  cm 

s 

17.0 

16.6 

15.3 

14.8 

12.2 

10.3 

U,  cm/sec 

5.54 

8.06 

9.45 

10.61 

12.95 

14.10 

-z  ’  /h 

s 

u/U 

0.00 

0.540 

0.540 

0.540 

0.540 

0.540 

0.540 

0.540 

0.02 

0.470 

0.440 

0.340 

0.420 

0.400 

0.420 

0.415 

0.04 

0.390 

0.200 

0.210 

0.270 

0.280 

0.370 

0.286 

0.06 

0.310 

0.120 

0.140 

0.160 

0.180 

0.320 

0.204 

0.10 

0.200 

0.070 

0.070 

0.100 

0.135 

0.235 

0.135 

0.15 

0.120 

0.035 

0.020 

0.050 

0.080 

0.170 

0.079 

0.20 

0.080 

0.020 

-0.010 

0.010 

0.040 

0.130 

0.045 

0.30 

0.000 

-0.010 

-0.035 

1.020 

-0.010 

0.070 

0.000 

0.40 

-0.045 

-0.025 

-0.040 

-0.020 

-0.050 

0.000 

-0.030 

0.50 

-0.075 

-0.040 

-0.040 

-0.020 

-0.080 

-0.060 

-0.052 

0.60 

-0.095 

-0.050 

-0.040 

-0.020 

-0.100 

-0.110 

-0.070 

0.70 

-0.105 

-0.055 

-0.040 

-0.020 

-0.120 

-0.140 

-0.080 

0.80 

-0.110 

-0.055 

-0.040 

-0.020 

-0.110 

-0.140 

-0.079 

0.85 

-0.100 

-0.050 

-0.040 

-0.020 

-0.100 

-0.130 

-0.063 

0.90 

-0.090 

-0.040 

-0.035 

-0.020 

-0.080 

-0.100 

-0.061 

0.95 

-0.060 

-0.035 

-0.020 

-0.010 

-0.050 

-0.070 

-0.041 

1.00 

-0.000 

-0.000 

-0.000 

-0.000 

-0.000 

-0.000 

-0.000 

Table  4 

Velocity  Distribution  in  Freshwater  Current 


tV  /o 

s 

0.02 

0.04 

u/U 

0.06 

Mean 

z 1  /h 

w 

0.00 

0.540 

0.540 

0.540 

0.540 

0.02 

0.668 

0.680 

0.676 

0.675 

0.04 

0.775 

0.771 

0.786 

0.777 

0.06 

0.853 

0.845 

0.865 

0.854 

0.08 

0.908 

0.882 

0.906 

9.898 

0.10 

0.948 

0.911 

0.937 

0.932 

0.12 

0.971 

0.935 

0.960 

0.955 

0.16 

1.012 

0.965 

0.995 

0.991 

0.20 

1.035 

0.989 

1.015 

1.013 

0.30 

1.055 

1.020 

1.040 

1.037 

0.40 

1.061 

1.035 

1.052 

1.049 

0.50 

1.061 

1.042 

1.052 

1.049 

0.60 

1.052 

1.042 

1.048 

1.047 

0.70 

1.032 

1.032 

1.032 

1.032 

0.80 

1.003 

1.028 

1.015 

1.015 

0.90 

0.970 

1.010 

0.995 

0.992 

1.00 

0.935 

0.991 

0.968 

0.965 

1.00 


Table  3 

Velocity  Distribution  In  Freshwater  Current 


Apg/p  = 

0.040;  h 

s 

+  h  *  30. 
w 

,  3  f.m 

Run 

18 

19 

20 

21 

22 

23 

Mean 

h  ,  cm 
w 

13.3 

16.6 

15.3 

14.8 

12.2 

10.3 

U,  cm/ sec 

z'/h 

w 

5.54 

8.06 

9.45 

10.61 

u/U 

12.95 

14.10 

0.00 

0.54 

0.54 

0.54 

0.54 

0.54 

0.54 

0.540 

0.02 

0.68 

0.67 

0.70 

0.72 

0.66 

0.65 

0.680 

0.04 

0.80 

0.79 

0.79 

0.81 

0.75 

0.73 

0.771 

0.06 

0.88 

0.86 

0.85 

0.85 

0.81 

0.81 

0.845 

0.08 

0.92 

0.91 

0.89 

0.89 

0.84 

0.84 

0.882 

0.10 

0.95 

0.94 

0.92 

0.91 

0.88 

0.86 

0.911 

0.12 

0.97 

0.95 

0.94 

0.93 

0.91 

0.90 

0.935 

0.16 

1.00 

0.98 

0.98 

0.96 

0.95 

0.92 

0.965 

0.20 

1.02 

0.99 

1.00 

0.98 

0.98 

0.96 

0.989 

0.30 

1.04 

1.02 

1.04 

1.01 

1.01 

1.01 

1.020 

0.40 

1.04 

1.03 

1.05 

1.03 

1.03 

]  .03 

1.035 

0.50 

1.04 

1.04 

1.06 

1.04 

1.04 

1.04 

1.042 

0.60 

1.03 

1.04 

1.06 

1.05 

1.04 

1.05 

1.042 

0.70 

1.01 

1.03 

1.05 

1.04 

1.04 

1.05 

1.032 

0.80 

1.00 

1.02 

1.04 

1.03 

1.04 

1.05 

1.028 

0.90 

0.97 

1.00 

1.02 

1.01 

1.02 

1.04 

1.010 

1.00 

0.95 

0.98 

1.00 

0.99 

1.00 

1.03 

0.991 

Table  2 

Dependence  of  Rise  of  Interface  and  Fall  of  Mater  Surface 
on  Densimetric  Froude  Number 


Run 

Ap 

10^  d  'j  /dx 

104U2/gH 

w 

104U2/gHQ 

104dh  /dx 
s 

102  x  u2/~ 

P 

34 

0.02 

3.75 

4.69 

2.21 

2.92 

2.34 

40 

7.36 

13.60 

5.96 

5.6 

7.80 

41 

10.95 

18.50 

8.76 

8.3 

9.25 

42 

20.91 

33.10 

16.98 

14.5 

11.55 

43 

30.00 

41.80 

24.08 

21.5 

20.70 

44 

33.40 

43.40 

28.6 

26.0 

21.70 

43 

42.40 

48.80 

36.2 

25.0 

24.40 

10 

0.04 

24.2 

34.6 

17.4 

10.0 

8.72 

1 1 

29.3 

47.0 

26.0 

15.0 

11.7 

12 

33.8 

47.0 

26.0 

15.0 

11.7 

13 

47.1 

61.8 

36.3 

18.0 

15.8 

14 

55.4 

67.8 

45.1 

24.5 

16.9 

16 

72.6 

83.1 

58.8 

30.5 

20.8 

17 

13.7 

20.5 

9.38 

8.0 

5.12 

24 

0.06 

11.8 

21.9 

9.80 

4.5 

3.84 

25 

23.4 

47.1 

21.6 

6.6 

7.85 

26 

43.1 

73.4 

34.9 

10.0 

12.28 

27 

55.4 

94.4 

47.7 

13.4 

15.70 

28 

72.0 

107.2 

61.2 

17.8 

17.80 

29 

82.4 

110.5 

71.4 

21.8 

18.40 

30 

92.2 

117.0 

82.1 

26.0 

19.50 

Table  1 

Rates  of  Rise  of  Interface  and  Fall  of  Water  Surface 


33 

O 

II 

1  & 
CO 

+  h  +  5 ;  H1  - 
w  1 

-  30. 

3  cm 

Run 

Ap/p 

g,  cm^/sec 

h  ,  cm 
s 

dh 

s 

,  4  4 

/dx  x  10 

dfi/dx  x  106 

39 

0.020 

37.  1 

17.8 

2.9 

3.75 

40 

55.7 

16.9 

5.6 

7.36 

41 

74.2 

15.8 

8.3 

10.45 

42 

111.3 

14.7 

14.5 

20.91 

43 

148.4 

12.8 

21.5 

30.00 

44 

185.6 

10.2 

26.0 

33.40 

45 

235.7 

7.7 

25.0 

42.40 

10 

0.040 

111.3 

14.9 

10.0 

24.2 

11 

148.4 

13.5 

15.0 

29.3 

12 

148.4 

13.5 

15.0 

33.8 

13 

185.6 

12.5 

18.0 

47.1 

14 

235.7 

10.0 

24.5 

55.4 

16 

285.7 

8.8 

30.8 

72.6 

17 

74.2 

16.3 

8.0 

13.7 

24 

0.060 

74.2 

16.6 

4.5 

11.8 

25 

111.3 

16.4 

6.6 

23.4 

26 

148.4 

15.8 

10.0 

43.1 

27 

185.6 

14.8 

13.4 

55.4 

28 

235.7 

12.9 

17.8 

72.0 

29 

285.8 

10.7 

21.8 

82.4 

30 

335.8 

8.9 

22.0 

92.2 

Note:  h  measured  at  station  x  =  2,000  . 
s 
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The  purpose  of  these  requirements  is  to  see  if  the  conditions  at  the  ends  have 

any  hearing  on  the  mixing  across  the  interface.  A  channel  of  30-m  length  and 

of  40-em  height  suffices.  To  minimize  the  uncertainties  relating  to  the  ver¬ 
tical  wall  stresses,  it  is  desirable  that  the  channel  be  wide,  say,  at  least 
90  cm.  1'lie  liquid  height  in  channel  may  be  chosen  as  30  cm,  with  an  underly¬ 
ing  saline  layer  depth  of  15  or  20  cm.  The  saline  water  may  be  introduced 

from  a  special  tank,  the  inflow  into  channel  being  either  from  the  bottom 
vertically  upward  or  from  the  channel  wall  laterally  at  the  level  of  the 
bottom,  the  flow  distribution  being  uniform  in  both  cases  along  the  entire 
length  of  the  bottom.  Care  is  taken  to  measure  accurately  the  quantity  of 
saline  inflow  to  assure  a  correct  balance  with  the  transport  of  saline  water 
across  the  interface.  The  densities  may  be  measured  by  parallel  wire  elec¬ 
trodes  in  one  of  the  following  manners.  The  liquids  from  different  levels 
may  be  sucked  and  collected  in  separate  bottles  in  the  manner  effected  bv 
Loftquist  (1960).  The  densities  of  the  collected  samples  in  bottles  can  be 
measured  in  later  times.  Or  the  electrode  may  be  attached  to  a  gage  which 
may  be  moved  vertically  through  the  entire  depth  for  the  local  measurement 
of  densities.  In  the  latter  case,  it  is  desirable  to  calibrate  the  electrode 
in  moving  saline  water.  This  is  a  matter  of  a  long  investigation.  The  liquid 
velocities  may  be  measured  by  deflecting  ribbons  as  in  the  present  tests.  To 
avoid  errors  of  uncertain  zero  readings,  the  deflections  may  be  determined 
through  strain  gage  indications,  the  gages  being  attached  to  the  ribbon  at  its 
clamped  area.  To  show  that  the  method  is  feasible  also  a  separate  investiga¬ 
tion  would  be  required.  The  fall  of  water  surface  may  be  measured  by  five 
optical  manometers  of  the  type  used  in  the  present  study.  For  accuracy,  how¬ 
ever,  it  is  necessary  that  each  manometer  is  provided  with  individual  tele¬ 
scopes.  1'hen,  the  front  cells  of  the  manometer  are  connected  to  a  tank  or  a 
trough  channel  of  constant  surface  level  attached  to  an  experimental  channel 
and  the  hinder  cells  to  the  channel  at  five  distant  points  with  equal  spacing 
between  them. 
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In  this  study  then  as  Figure  38  shows,  the  flow  adjacent  to  the  bottom  is 
laminar  and  the  stress  would  be 


T 

O 


u 


m 


or  in  terms  of  h  and  U 
s 


t  2u  h  / Uh  Y 

o  _  _ m  __s  /  s  l 

„2  U  d  l  v  / 

pU  m  \  / 


Using  the  values  of  u  /U  and  h  /d  mentioned  above 

m  s  m 


(117) 


(118) 


Comparing  with  Equation  112  it  is  seen  that  the  estimate  of  the  stress  from  the 
observed  velocities  is  ten  times  smaller  than  that  computed  from  the  hydro- 
dynamic  Equation  101. 

56.  The  severe  disagreement  between  the  two  determinations  of  the 
bottom  stress  above  indicated  is  disappointing.  The  cause  of  the  disagree¬ 
ment  in  one  part  is  the  difficulty  of  low  flow  measurements  and  in  another 
part,  the  inaccuracies  in  the  measurements  of  layer  depths  and  slopes  of 
interface  and  of  upper  free  surface.  Also  the  method  of  evaluating  the 
vertical  wall  stresses  is  open  to  question.  Further,  probably  Equation  102 
does  not  take  into  account  the  actual  saline  flow  conditions  met  with  in  the 
tests.  It  will  be  remembered  that  to  balance  the  saline  water  transfer  across 
the  interface,  saline  flows  were  introduced  laterally  at  the  channel's  two 
extreme  ends.  The  dynamic  effect  of  these  flows  on  the  saline  layer  was  not 
considered  in  developing  Equation  102. 


Recommendations  for  a  New  Research 


57.  In  view  of  the  fact  that  the  dependence  of  the  transport  of  saline 
water  across  the  interface  upon  densimetric  Froude  number  is  not  in  agree¬ 
ment  with  the  results  of  other  investigations,  it  is  desirable  that  the 
problem  be  considered  in  a  modified  channel  specially  in  order  to  pay  atten¬ 
tion  to  a  possible  mixing  in  the  area  of  the  freshwater  entrance.  In  the  new 
channel,  the  bottom  is  inclined  at  the  entrance  and  the  efflux  is  over  a  weir 
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Figure  S8.  Dependence  of  bottom  stress  on  Reynolds  number 
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in  terms  of 


and  li  the  criterion  may  now  be  written  also  as 
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900 
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LI 


d  u 
m  m 


(1 


From  Figure  1.5,  representing  the  velocities  of  the  salt  layer,  one  may  take 

u  /L'  =  0.15  and  d  /h  =  0.3  .  I'sing  these  values,  the  condition  for  the 
m  ms 

luminary  regime  in  the  part  adjacent  to  the  bottom  is 


Uh 
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-  41,500 
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0.066 


(110) 
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which  is  equally  useful  as  the  logarithmic  form  and  is  based  on  scattered 
data  . 

54.  For  the  evaluation  of  the  interfacial  stress  when  the  fall  of 
the  water  surface  is  not  obtainable,  one  may  proceed  supposing  that  the 
momentum  effects  in  the  saline  water  are  negligible.  Pedersen  proceeded  on 
this  basis  to  obtain  data  from  various  authors  to  establish  the  relation 
presented  by  Equation  107.  To  proceed  in  the  same  manner  for  the  present 
investigation,  one  obtains  by  eliminating  d(S/dx  from  Equations  92  and  98 
after  ignoring  the  momentum  effects  in  the  second  equation,  the  result 


gh  Ap 
w  s 


dh  h 

— s  -  1  -w 
dv  o  B 


1  +  h  , 
w  s 


dll) 


Evaluations  made  on  this  basis  are  shown  in  the  fourth  column  of  Table  25. 

The  average  interfacial  shear  stress  thus  evaluated  is  nearly  30  percent  higher 
than  when  the  bottom  stress  and  momentum  effects  are  not  ignored. 

55.  In  Figure  38  the  data  of  the  bottom  stress,  which  are  given  in  the 
fifth  column  of  Table  25,  are  plotted  against  the  Reynolds  number  Uh  /v  ■ 
Although  the  scatter  of  the  data  points  is  severe,  nevertheless  the  average 
distribution  may  be  represented  by  the  straight  line  drawn  with  the  equation 
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(112) 


The  relevance  of  this  result  with  the  measured  velocities  in  the  vicinity  of 
the  bottom  may  be  next  examined.  Relying  on  a  result  from  open  channel  flow, 
the  criterion  for  the  laminar  regime  in  the  layer  adjacent  is 

_  d 

u  —  <  500  (113) 

\) 


where  u  is  the  average  velocity  in  the  layer,  and  d^  the  distance  of  the 

maximum  u  from  the  bottom.  For  a  laminar  layer 
m 


67 


( 

I 


Figure  37.  Dependence  of  interfacial  stress  on  Reynolds  number 

overflow.  Aiming  to  obtain  a  relation  applying  solely  to  an  overflow,  we  place 

U  =  1.79  (U  -  U.)  and  identify  y  -  y  with  the  freshwater  depth  h  and 
m  t  1  o  Kw 

next  writing 

X./2  =  t  /oU2  (108) 

1  s 

R  =  Uh  /v 
w 

the  Pedersen  relation  becomes 


3.45  [ln(RVyi)-  lo] 


(109) 


The  curve  shown  in  Figure  37  is  from  this  relation,  and  the  data  points  of 
the  present  investigation  are  in  fair  agreement  with  the  Pedersen  formulation. 
The  equivalent  formula  in  a  power  form  is 
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Table  10 

Interface  Layer  Effective  Thickness,  1 


(Observed  at  x  =  3,000  cm) 


0.02 


0.04 


0.06 


Run 

U,  cm/sec 

h  ,  cm 
w 

1  ,  cm 

h  /I 
w 

39A 

2 . 64 

14.1 

0.95 

14.8 

40A 

4.00 

13.9 

0.90 

15.5 

4  1 A 

5.60 

13.3 

0.85 

15.6 

^2A 

8.12 

13.8 

1.00 

13.8 

4  3A 

10.58 

14.0 

1.40 

10.0 

44  A 

11.72 

15.8 

0.90 

17.5 

4  5  A 

13.60 

17.3 

1.40 

12.2 

]. 

2.97 

12.5 

1.15 

10.8 

2 

3.42 

11.3 

1.30 

8.8 

3 

6.58 

11.3 

0.91 

12.4 

4 

7.15 

15.6 

1 .08 

12.5 

5 

10.30 

14.4 

1. 10 

13.1 

6 

12.70 

14.6 

1.16 

12.6 

7 

15.40 

15.3 

1.50 

10.1 

8 

17.85 

16.1 

0.82 

19.7 

24 

5.85 

12.9 

0.75 

17.2 

25 

8.65 

12.9 

0.65 

19.9 

26 

11.70 

12.7 

0.70 

18.1 

27 

14.30 

13.0 

0.80 

16.4 

28 

16.52 

14.3 

0.90 

15.9 

29 

17.70 

16.0 

1.10 

14.4 

30 

19.80 

17.0 

1.35 

12.5 

f 


Table  11 

Interface  Layer  Effective  Thickness,  1 
(Observed  at  x  =  2,000  cm) 


■V>„/p 

3 

Run 

U,  cm/ sec 

h  ,  cm 
w 

1  ,  cm 
u 

h  / 1 
w  u 

0.02 

39B 

2.58 

14.6 

2.0 

7.30 

40B 

3.90 

14.3 

2.3 

6.18 

41B 

5.12 

14.4 

2.1 

7.42 

42B 

7.31 

15.3 

2.9 

5.80 

43B 

9.28 

16.0 

3.1 

5.10 

44B 

9.42 

19.3 

3.3 

5.80 

45B 

10.95 

21.5 

3.6 

6.02 

0.04 

9 

10.38 

17.2 

2.2 

7.70 

18 

5.54 

13.4 

2.0 

6.70 

19 

8.06 

13.8 

1.9 

7.18 

20 

9.45 

15.7 

2.5 

6.28 

21 

10.67 

17.4 

2.4 

7.30 

22 

12.95 

18.2 

2.8 

6.58 

23 

14.10 

20.3 

3.6 

5.72 

0.06 

24A 

5.36 

14.0 

1.8 

7.84 

25A 

7.96 

14.0 

2.2 

6.30 

26A 

10.24 

16.3 

2.5 

6.54 

27A 

12.30 

15.1 

2.5 

6.02 

28A 

14.28 

16.5 

3.1 

5.26 

29A 

15.78 

18.1 

3.0 

6.00 

30A 

16.70 

20.1 

3.5 

5.60 

Table  12 


Density  and  Salinity 

Relations 

for  Saltwater  Solutions 

p  •=  Density 

of  Water; 

p  +  Ap  =  D 

ensity  of 

Saline  Solution 

p/p  =  0.465 
m 

Ap/p 

Ut 

S 

€ 

S/(Ap/p 

0.0053 

1.89 

0.01 

1.032 

1.89 

0.0125 

1.62 

0.02 

0.826 

1.60 

0.0268 

1.54 

0.04 

0.756 

1.49 

0.0413 

1.51 

0.06 

0.734 

1.45 

0.0559 

1.51 

0.08 

0.734 

1.43 

0.0707 

1.52 

0.10 

0.734 

1.41 

0.0857 

1.52 

0.12 

0.744 

1.40 

0.1009 

1.52 

0.14 

0.748 

1.40 

Table  13 

Saline  Density  Distribution  in  Fresh  Water 

2 

Run  27B,  Apg/Ap  =  0.06;  q  =  185.6  cra~/sec 
x  =  1,000  cm;  U  =  11.4  cm/sec 


z ,  ft 

Ap/Ap 

z,  ft 

Ap/Ap 

-0.095 

0.970 

0.010 

0.159 

-0.055 

0.970 

0.014 

0.00734 

-0.035 

0.970 

0.015 

0.00648 

-0.025 

0.941 

0.019 

0.00546 

-0.016 

0.922 

0.0024 

0.00444 

-0.006 

0.836 

0.034 

0.00375 

-0.005 

0.789 

0.044 

0.00341 

-0.001 

0.624 

0.064 

0.00239 

0.000 

0.508 

0.104 

0.00188 

0.004 

0.193 

0.204 

0.00102 

0.005 

0.162 

0.304 

0.00086 

Table  14 

Saline  Density  Distribution 

in  Fresh  Water 

Run  27B,  Ap/Ap 

s 

=  0.06;  q  = 

2 

185.6  cm  / sec 

x  =  3,000 

cm;  U  =  14 

. 3  cm/ sec 

z ,  cm 

Ap/Ap 

s 

z ,  cm 

Ap/Aps 

-0.088 

0.895 

0.037 

0.0123 

-0.048 

0.886 

0.052 

0.0097 

-0.028 

0.861 

0.062 

0.00846 

-0.018 

0.844 

0.072 

0.00665 

-0.007 

0.768 

0.092 

0.00511 

-0.003 

0.613 

0.112 

0.0046 

0.002 

0.448 

0.162 

".00358 

0.007 

0.261 

0.212 

0.00239 

o.o:_2 

0.0524 

0.262 

0.00175 

0.017 

0.0283 

0.312 

0.00136 

0.022 

0.0165 

0.412 

0.00136 

0.030 

0.0146 

Table  15 

Saline  Density  Distribution  In  Fresh  Water 
Run  27B,  Ap^/Ap  =  0.06;  q  =  185.6  cm^/sec 

x  =  2,000  cm;  U  =  13.0  cm/sec 


z ,  ft 

Ap/APg 

-0.043 

0.959 

-0.033 

0.951 

-0.023 

0.933 

-0.013 

0.878 

-0.008 

0.848 

-0.003 

0.675 

0.002 

0.364 

0.007 

0.1039 

0.012 

0.0320 

0.015 

0.0156 

0.022 

0.0114 

z ,  ft 

Ap/Apg 

0.027 

0.0104 

0.032 

0.00906 

0.037 

0.00808 

0.047 

0.00586 

0.057 

0.00515 

0.077 

0.00319 

0.097 

0.00279 

0.147 

0.00195 

0.243 

0.00142 

0.347 

0.00080 

Table  16 

Saline  Transport  Across  Interface 

Run  44C,  Apg/p  =  0.02;  Run  9,  Ap  /p  =  0.04;  Run  27B,  Apg/p  =  0.06 


q  =  185.6  cm  /sec 


Ap  /p 

s 


10  L,  cm 


U ,  cm/ sec 


h  ,  cm 
w 


2Qs2.  cm  /sec 


0.02  1 
2 
3 


7.8 

21.3 

0.66 

8.7 

19.0 

1.06 

9.5 

15.3 

2.36 

0.04  1 

2 
3 


7.8 

19.3 

0. 

.59 

8.5 

17.4 

1. 

.13 

9.8 

14.5 

2. 

,15 

1  10.6  16.2 

2  11.4  14.2 

3  12.0  13.2 


0.36 

0.79 

1.59 


0.06 


Table  17a 


Ap  ,/p  =  0.02;  L  =  2,000  cm;  v  =  0.0091 
U  =  q/h 

w 


Run 

h  ,  cm 
w 

U,  cm/sec 

10  u  _ 

ml 

10  u  „ 

mz 

39B 

15.8 

2.35 

3.36 

0.06 

40B 

15.7 

3.36 

4.52 

0.18 

4  IB 

15.5 

4.78 

5.72 

0.25 

42B 

16.3 

6.86 

9.40 

1.00 

43B 

17.5 

8.52 

13.62 

1.90 

44B 

20.6 

9.05 

13.70 

2.48 

44C 

20.7 

9.04 

12.60 

2.90 

45B 

22.4 

10.50 

18.90 

4 . 66 

Table  17b 


Ap  /p  =  0.02;  L  =  3,000  cm;  v  =  0.0091 
s _ 

U  =  q/h 

w 


Run 

h  ,  cm 
w 

U,  cm/ sec 

10  U  . 
ml 

10  U  9 

m2 

39A 

14.6 

2.54 

3.31 

0.03 

40A 

14.5 

3.84 

6.00 

0.13 

41A 

14 . 6 

5.05 

7.10 

0.48 

42A 

16.0 

7.05 

10.4 

1.48 

43A 

17.8 

8.35 

19.0 

2.58 

44A 

20.6 

9.00 

19.8 

3.60 

44C 

20.6 

9.00 

16.2 

3.63 

45A 

23.0 

10.00 

21.4 

3.50 

Table  18a 


Apg/p  =  0.04;  L  =  2,000  cm;  v  =  0.0108 
U  =  q/h 

w 


Run 

h  ,  cm 
w 

V ,  cm/ sec 

18 

14.5 

5.08 

19 

14.9 

7.40 

20 

16.4 

8.68 

21 

19.2 

9.70 

22 

19.7 

11.90 

23 

22.0 

12.95 

4  4 

10  1)  ,  ,  cm/sec  10  U  „  ,  cm/sec 

nu  m2 


7.90 

0.09 

12.60 

0.05 

15.60 

0.98 

13.80 

2.83 

15.30 

4.53 

22.80 

8.60 

Table  18b 


Apg/p  =  0.04;  L  =  3,000  cm;  v  =  0.0106 
U  =  q/hw 


Run 

h  ,  cm 
w 

U,  cm/sec 

4 

10  U  _ ,  cm/ sec 
ml 

4 

10  U  cm/ sec 

m/ 

1 

14.2 

1.88 

3.19 

0.01 

2 

13.6 

3.06 

5.60 

0.08 

3 

13.9 

4.16 

7.60 

0.28 

4 

17.8 

6.26 

10.80 

0.29 

5 

17.8 

8.36 

14.70 

1.68 

6 

18.0 

10.32 

16.50 

3.70 

7 

19.6 

12.95 

25.40 

5.80 

8 

20.3 

14.15 

27.80 

8.18 

Table  19a 


Aps/p  =  0.06;  L  =  2,000  cm;  v  =  0.0098 


Run 

h  ,  cm 
w 

U,  cm/ sec 

10  U  _ ,  cm/sec 
ml 

10  U  cm 

mz 

24  A 

L  5 . 2 

4.85 

7.41 

0.006 

2  5A 

15.2 

7.30 

10.00 

0.017 

26A 

15.5 

9.40 

13.52 

0.079 

27A 

16.4 

11.30 

19.80 

2.39 

2  8  A 

17.8 

13.10 

19.00 

3.45 

29A 

1.9.7 

14.50 

23.10 

5.16 

30  A 

21.8 

15.35 

30.40 

7.90 

Table  19b 

Ap  / p  =  0.06; 

L  =  3,000  cm;  v  =  0.0103 

q  =  Uh 

w 

Run 

h  ,  cm 
w 

U,  cm/sec 

10^U  - ,  cm/sec 
ml 

10^U  ~ ,  cm/sec 
mz 

24 

14.15 

5.06 

7.9 

0.088 

25 

14.4 

7.62 

11.4 

0.37 

26 

14.8 

9.78 

15.5 

1.75 

27 

15.9 

11.75 

19.6 

3.85 

28 

17.8 

13.65 

21.4 

5.12 

29 

19.5 

15.10 

23.30 

7.10 

30 

21.0 

15.95 

29.80 

10.10 

Table  20a 


Ap  /p  =  0.02;  L  =  2,000  cm;  -  =  0.0091 


Run 

104  ^ 

“>4¥ 

39B 

1.42 

0.02 

40B 

1.35 

0.05 

41B 

1.12 

0.05 

42B 

1.37 

0.14 

43B 

1.60 

0.22 

44B 

1.44 

0.27 

44C 

1.39 

0.32 

45B 

1.81 

0.44 

104  j~ 

Um 

2  2 

10  F, 

A 

io3r 

1.44 

1.77 

4.07 

1.40 

3.69 

6.10 

1.17 

7.40 

8.10 

1.51 

14.82 

12.20 

1.82 

21.16 

16.28 

1.71 

20.25 

20.30 

1.71 

20.25 

20.30 

2.25 

25.10 

25.80 

Table  20b 


Apg/p  =  0.02;  L  =  3,000  cm;  v  =  0.0091 


Run 

104 

io4  -2^- 
U 

39A 

1.33 

0.01 

40A 

1.65 

0.03 

41A 

1.37 

0.09 

42A 

1.35 

0.21 

43A 

1.81 

0.31 

44A 

1.66 

0.40 

44C 

1.73 

0.40 

45A 

1.61 

0.35 

104  !T 

Um 

2  2 

10  fa 

3 

10  R 

1.34 

2.02 

4.07 

1.68 

4.84 

6.14 

1.46 

8.53 

8.04 

1.56 

15.8 

12.30 

2.12 

19.9 

16.20 

2.06 

19.8 

20.80 

2.13 

19.8 

20.80 

1.96 

24.1 

26.00 

Table  21a 


Apg/j  =  0.04;  L  -  2,000  ctn;  v  =  0.0108 


Run 

,  u 

u 

io4 

18 

1  .  5b 

0.02 

19 

1.  ,0 

0.01 

20 

L .  80 

0.11 

21 

13.8 

0.29 

22 

1.28 

0.38 

23 

1.76 

0.43 

2  2 

10  F, 

A 

103R 

1.58 

4.54 

6.08 

1.71 

9.36 

10.32 

1.91 

11.7 

13.75 

1.67 

12.4 

17.22 

1.66 

18.2 

17.22 

2.19 

19.3 

21.80 

Table  21b 

Apg/p  =  0.04;  L  =  3,000  cm;  v  =  0.0106 


Run 

/  U  1 

10  Hr 

1 

1.69 

0.01 

2 

1.82 

0.03 

3 

1.82 

0.06 

4 

1.71 

0.05 

3 

1.75 

0.21 

6 

1.59 

0.36 

7 

1.95 

0.44 

8 

1.95 

0.58 

i»‘r 

2  2 

10  F* 

103R 

1.70 

0.64 

3.50 

1.85 

1.99 

5.28 

1.88 

3.13 

7.02 

1.76 

5.66 

10.53 

1.96 

9.98 

14.02 

1.95 

15.2 

17.57 

2.35 

18.9 

22.30 

2.53 

25.0 

27.10 

Table  22a 


Apg/p  =  0.06;  L  =  2,000  cm;  v  =  0.0098 


Run 

io4  Y 

io4  ^ 
u 

24A 

1.53 

0.00 

25A 

1.36 

0.00 

26A 

1.44 

0.01 

27A 

1.95 

0.21 

28A 

1.45 

0.26 

29A 

1.59 

0.37 

30A 

1.98 

0.62 

i»4r 

2  2 

10  F, 

A 

3 

10  R 

1.53 

2.66 

7.56 

1.36 

6.15 

11.34 

1.45 

9.67 

15.10 

1.96 

14.7 

18.92 

1.71 

16.5 

24.00 

1.96 

18.1 

29.10 

2.50 

18.4 

34.28 

Table  22b 

Apg/p  =  0.06;  L  =  3,000  cm;  v  =  0.0103 


Run 

“4-f- 

i°4-r 

61 

o 

r— 1 

2  2 
io  f: 

A 

3 

10  R 

24 

1.55 

0.02 

1.57 

2.99 

7.30 

25 

1.50 

0.05 

1.55 

6.86 

10.95 

26 

1.60 

0.18 

1.78 

10.24 

14.10 

27 

1.67 

0.33 

2.01 

14.1 

18.25 

28 

1.56 

0.38 

1.94 

17.6 

23.20 

29 

1.59 

0.52 

2.11 

19.9 

28.00 

30 

1.87 

0.63 

2.50 

19.8 

33.00 

Table  23 

Measured  Quantities  Relating  to  the  Stresses 


h  =  H  -  h  ;  H  =  30.3  cm;  h  Measured  at  x  =  2,000  cm 
w  s  s 


Run 

h  ,  cm 
s 

h  ,  cm 
w 

U,  cm/sec 

dh  / dx 
s 

d5/dx 

39 

0.02 

17.8 

12.5 

2.97 

2.90  x  10"4 

3.75  x  io-6 

90 

16.9 

13.4 

4.22 

5.60 

7.36 

91 

15.8 

14.5 

5.12 

8.30 

10.9 

92 

14.7 

15.6 

7.12 

14.5 

20.9 

93 

12.8 

17.5 

8.48 

21.5 

30.0 

49 

10.2 

20.1 

9.29 

26.0 

33.4 

45 

7.7 

22.6 

10.42 

35.0 

42.4 

n 

0.04 

14.9 

15.4 

7.22 

10.0  x  io-4 

24.2  x  lo~6 

1 1 

13.3 

16.8 

8.90 

15.0 

29.2 

12 

13.5 

16.8 

8.80 

15.0 

35.6 

1  J 

12.5 

17.8 

10.42 

18.0 

47.1 

14 

10.0 

20.3 

11.60 

24.5 

55.4 

lb 

8.8 

21.5 

13.26 

3.05 

72.6 

1  7 

16.3 

14.0 

5.30 

8.0 

13.7 

24 

0.06 

16.6 

13.7 

5.42 

4.5  x  io"4 

11.8  x  io"6 

25 

16.4 

13.9 

8.00 

6.6 

23.4 

26 

15.8 

14.5 

10.23 

10.0 

43.1 

27 

14.8 

15.5 

11.95 

13.4 

55.4 

28 

12.9 

17.4 

13.55 

17.8 

72.0 

29 

10.7 

19.6 

14.58 

21.8 

82.4 

30 

8.9 

21.4 

15.67 

21.0 

92.2 

Table  24 


Coefficient  of  Resistance  of  the  Channel  Sidewalls 
2b  =  22.7  cm 


(U  and 

in  cgs  units) 

Run 

.vs/,, 

103v 

1,1/4 

3 

10  X 

o 

103\  h  /B 
o  w 

103  dh  , 
s 

39 

0.020 

9.39 

1.31 

5.55 

3.40 

0.29 

40 

1.43 

5.10 

3.16 

0.56 

41 

1.50 

4.86 

3.10 

0.83 

42 

1.54 

4.73 

3.35 

1.45 

43 

1.71 

4.27 

3.29 

2.25 

44 

1.73 

4.20 

3.70 

2.60 

45 

1.79 

4.08 

4.08 

3.50 

10 

0.040 

11.11 

1.64 

4.64 

3.09 

1.00 

11 

1.72 

4.41 

3.18 

1.50 

12 

1.72 

4.41 

3.18 

1.50 

13 

1.79 

4.25 

3.50 

1.80 

14 

1.85 

4.11 

3.63 

2.45 

16 

1.91 

3.98 

3.88 

3.05 

17 

1.51 

6.00 

3.60 

0.80 

24 

0.060 

10.07 

1.53 

4.89 

2.98 

0.45 

25 

1.68 

4.46 

2.72 

0.66 

26 

1.79 

4.19 

2.70 

1.00 

27 

1.86 

4.03 

2.78 

1.34 

28 

1.91 

3.93 

3.06 

1.78 

29 

1.95 

3.84 

3.29 

2.18 

30 

1.99 

3.76 

3.50 

2.60 

Table  25 


I 


Experimental  Evaluation  of  Interfaci al  and  Bottom  Stresses 
F.  =  17  yj  (2-  s/.  )ghs 


Run 

/. 

s 

"r  / 

s 

;  U" 

7  / 

S 

r  /.  u2 

0 

U  h  / 
w 

1  OF , 

39 

0.02 

1.48 

••  10~3 

2 . 29 

'  io-3 

3.32  x  10'3 

3.93  ••  103 

1.53 

40 

1.70 

2.00 

2.20 

5.90 

2.79 

41 

1.99 

2.61 

1.90 

7.90 

3.03 

42 

1.64 

1.92 

1.11 

11.92 

3.40 

43 

1.74 

2.12 

1.19 

14.95 

4.55 

44 

1.4 

1.92 

0.75 

19.85 

4.66 

45 

1.16 

2.30 

0.67 

25.10 

4.94 

10 

0.04 

2.99 

•  io“3 

3.86 

*  10“ 3 

2.80  x  10-3 

9.78  x  103 

2.45 

11 

1.54 

3.56 

2.42 

13.02 

3.42 

12 

2.50 

3.56 

2.42 

13.02 

3.42 

13 

1.34 

1.81 

1.19 

15.25 

3.98 

14 

2.10 

2.77 

0.80 

25.10 

4.11 

16 

1.79 

2.25 

1.09 

25.90 

4.56 

17 

2.48 

6.42 

6.64 

2.26 

24 

0.05 

2.08 

X  10"3 

4.98 

*  10-3 

2.64  x  10~3 

m 

o 

r—4 

o 

rv. 

1.96 

25 

1.58 

2.75 

2.08 

10.00 

2.90 

2b 

2.16 

2.31 

1.49 

14.50 

3.51 

27 

1.74 

2.10 

1.11 

18.00 

3.96 

28 

1.88 

2.16 

1.01 

23.00 

4.22 

29 

1.95 

3.35 

0.81 

28.60 

4.29 

30 

1.81 

2.11 

0.61 

32.50 

4.41 

i 


APPENDIX  A:  NOTATION 


Wfdt  of  ribbon;  or  half  channel  width 
Width  of  channel,  B  =  2b 

A  constant  relating  to  interlace  slope.  Equation  L 

Drag  coefficient  of  ribbons 

Modulus  of  elasticity  of  ribbon  material 

O 

2  IT 

Dens  i  me  trie  Froucte  number:  F  =  — >- - . — r- - 

A  g (Ap  / p  )  h 
s  w 

Acceleration  due  to  gravity 

Depth  of  saline  water  layer;  interface  heights 
Depth  of  freshwater  layer 

Elevation  of  free  surface  measured  from  the  channel  bottom 

H  -  h  +  h 
s  w 

Elevation  of  free  surface  in  the  entrance  and  exit  segments  of 
clianne  1 

Optical  manometer  telescope  indication,  k  =  88  i  ^ 

Effective  thickness  of  interface  layer  based  on  densities  as 
defined  by  Equation  26 

Effective  thickness  of  interface  1  aver  based  on  velocities  as 
defined  hv  Equation  '3  5 

Length  of  interface  line  between  x  =  0  and  x  -  x;  length  of  wires 
IT  ossure 

langent ial  arcus  on  a  plane  normal  to  x-axts  ami  in  t he 
d  i  red  ion  el  ■.  -  i  i  s 
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! .  it.. lent  t  i '  •  i  e  ,  .a  i  plane  normal  to  r  axis  and  in  tlu  direct  ion 
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Kevnolds  number  I’ll  i  .  ;  or  telescope  scale  reaching;  or  cuwntun 

w 

of  velocity  meter  ribbon  at  distance  s  ;  or  effective  resist;!  toe 
tress  Richardson  number  defined  nv  Equation  li 
Interfaeial  Richardson  number  del  i  nod  bv  Equation  lb 
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s  Distance  measured  along  ribbon  length;  or  the  spacing  between  the 
wires  of  salinity  measuring  device;  or  the  distance  between  the 
cylinders  of  the  optical  manometer 

S  Distance  of  telescope  scale  from  the  manometer  cells;  or  salinity 
t  Thickness  of  ribbon 

u  Depth  averaged  of  u  for  a  given  y  along  the  depth  h  of  fresh 
water  layer 

u^w^  Horizontal  and  vertical  velocity  components,  respectively,  of  inter 
face  points 

u  Saline  water  mixing  velocity  at  the  interface 
m 

ug  Current  velocity  at  the  freshwater  free  surface 

u,w  Horizontal  and  vertical  velocities  at  a  point 

U  Average  value  of  u  ,  averaged  over  the  channel  width  B;  average 
velocity  of  flowing  fresh  water 

U  Threshold  critical  velocity  associated  with  the  initiation  of 
c  ,  . 

mixing 

U.,W,  Average  of  values  of  u.,w.  ,  respectively,  averaged  over  the 

1  channel  width  B  11 


ml 


U 


m2 

U 


Average  of  u  ,  averaged  across  the  channel  width  B  ; 

U  =  U  .  +  U  “ 
m  ml  m2 

Component  part  of  transport  defined  by  Q^/L  where  Qg^  denotes 

the  volume  of  salt  carried  by  the  interface  per  unit  time  and  L, 
the  interface  length 

Freshwater  component  defined  by  Q^/L  »  where  Qg2  denotes  the 
volume  of  salt  carried  by  fresh  water  per  unit  time 
Average  value  of  ug  ,  averaged  over  the  channel  width  B 
Horizontal  coordinate  drawn  in  the  direction  of  flow  of  fresh  water 


y  Horizontal  coordinate  drawn  laterally  with  respect  to  channel  walls 

z  Elevation  of  points  measured  from  the  channel  bottom 

z'  Elevation  of  points  in  the  freshwater  layer  measured  from  the 

interface  line  where  Ap/pg  =  0.5 

i.  Length  of  ribbon  of  velocity  measuring  device 

a  Boussinesq  coefficient  of  velocity  distribution  as  in  Equation  86; 
or  a  coefficient  relating  to  salinity  as  in  Equation  48;  also 
specific  ratio  by  mass  for  solutions  as  in  Equation  46 

8  Coefficient  of  velocity  distribution  in  saline  water  as  defined  by 
Equation  95;  or  a  coefficient  relating  to  density  distribution 
in  freshwater  layer  as  in  Equation  39 

6  Fall  of  freshwater  free  surface;  H  =  H  -  6 

o 

6  Relative  displacement  of  liquid  surfaces  in  the  two  cells  of  the 
optical  manometer 


A2 


5  Lower  bound  of  interfacial  layer  where  Ap/Ap  =  0.99 
o  s 

Ap'  Average  value  of  Ap  in  fresh  water,  averaged  over  depth  h^ 

Ap  Excess  of  density  of  saline  water  of  lower  layer  over  that  of 
fresh  water 

e  Absorption  contents  in  mixtures.  Equation  47 

A^  Friction  coefficient  of  interface  as  defined  by  Equation  108 

A^  Friction  coefficient  of  wall  of  part  in  contact  with  saline  water 

as  defined  by  Equation  99 

A  Friction  coefficient  of  wall  of  part  in  contact  with  fresh  water  as 

J  defined  by  Equation  92 

V  Coefficient  of  viscosity 

v  Kinematic  viscosity  of  water 

p  Density  of  fresh  water 

p  Density  of  solute 

m 

t  Specific  resistance  of  electrolyte  in  ohms.  Equation  36 

t  Bottom  frictional  stress;  r  average  value  averaged  over  the 
channel  width  B 

x  Tangential  stress  at  the  interface;  rg  average  value  averaged 
over  the  channel  width  B 

Frictional  stress  of  vertical  walls;  averaged  value;  averaged 

separately  over  parts  of  wall  in  contact  either  with  fresh  water 


or  saline 

water 

n 

Volume 

of 

water 

n 

m 

Volume 

of 

solute 

Volume 

of 

solution 
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